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ABSTRACT 

We present an analysis of the galaxy-scale gaseous outflows from the FIRE (Feedback in 
Realistic Environments) simulations. This suite of hydrodynamic cosmological zoom simu¬ 
lations resolves formation of star-forming giant molecular clouds to z = 0, and features an 
explicit stellar feedback model on small scales. Our simulations reveal that high redshift 
galaxies undergo bursts of star formation followed by powerful gusts of galactic outflows 
that eject much of the ISM and temporarily suppress star formation. At low redshift, however, 
sufficiently massive galaxies corresponding to L*-progenitors develop stable disks and switch 
into a continuous and quiescent mode of star formation that does not drive outflows far into 
the halo. Mass-loading factors for winds in L*-progenitors are rj ss 10 at high redshift, but 
decrease to rj <C 1 at low redshift. Although lower values of rj are expected as halos grow 
in mass over time, we show that the strong suppression of outflows with decreasing redshift 
cannot be explained by mass evolution alone. Circumgalactic outflow velocities are variable 
and broadly distributed, but typically range between one and three times the circular velocity 
of the halo. Much of the ejected material builds a reservoir of enriched gas within the cir¬ 
cumgalactic medium, some of which could be later recycled to fuel further star formation. 
However, a fraction of the gas that leaves the virial radius through galactic winds is never 
regained, causing most halos with mass M\ < 10 I2 Mq to be deficient in baryons compared to 
the cosmic mean by z = 0. 

Key words: galaxies: formation — galaxies: evolution — stars: formation — cosmology: 
theory 


1 INTRODUCTION 

Young, hot stars release energy and momentum in the form of stel¬ 
lar winds, ionizing photons, radiation pressure, and powerful super¬ 
nova explosions. Nearby gas absorbs this energy and momentum, 
is heated up and pushed away supersonically from the dense re¬ 
gions where stars formed, driving shocks that propagate through 
the interstellar medium (ISM). Working in unison, a population of 
young stars is capable of accelerating a substantial mass of gas be¬ 
yond the escape velocity of the local gravitational potential, totally 
disrupting the local star-forming cloud, and potentially driving a 
large-scale galactic wind ([ Mathews & Baker| 1971[ |~Larson ||1974[ 
Haehnelt 19951. The effects of these energetic processes, both on 

the local state of the star-forming gaseous clouds, and on the dy¬ 
namics and morphological structure of the galaxy as a whole, are 
collectively known as stellar feedback. 

Winds attributed to stellar feedback have been directly ob¬ 
served in galaxies near and far using a variety of techniques in every 
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window of the electromagnetic spectrum (see |Veilleux et al.|2005| 
and references therein for a review). The approaches can gener¬ 
ally be divided into either studying diffuse, heated gas following an 
outflow (e.g. |Heckman et al.|1990| l or direct detection of kinemat¬ 
ically distinct material from the rest frame of the host galaxy via 


Doppler shifted emission (e.g. MartinJ1998l and absorption lines 


(e.g. 

Heckman et al. 2000[ Grimes et al. 2009, Weiner et al. 2009, 

Steic 

el et al.|20101 Martin et al. 2012; Bouche et al.i|2012; Rubin 


|et al.|20l~4j . Taken together, observations reveal that while galactic 
winds are nearly ubiquitous around starbursting galaxies, the rela¬ 
tive strengths of the winds as measured by the amount of material 
they move can vary substantially from case to case. In addition the 
winds’ kinematic properties (speeds range from tens to thousand 
of km/s) and phase compositions (cold molecular gas, warm gas 
detectable in the UV, and hot X-ray gas can often be seen simulta¬ 
neously) are very diverse. 

The net result of these outflows is the transport of material 
from the dense star-forming regions to the so-called circumgalac¬ 
tic medium (CGM), which is loosely defined as the gas outside the 
galaxy but still within the gravitational potential of the dark matter 
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halo that the galaxy inhabits. While it has been difficult to directly 
probe the total mass and phase composition of the CGM, detections 
of metal absorption everywhere from a few kpc from galactic cen¬ 
ters to regions at their virial radius and beyond hints at a vast reser¬ 
voir of gas that has been at least partially enriched with metals from 
galactic outflows JChMte£jdJ[1998 2001jJGauthierrta2|2009J 


Tu^inson^et^lj[20n^ Kacprzak & Churchill||2011| Werk et al.| 

2014| Lehner et al.|2015|. The high redshift CGM has also been 

extensively observed ( [Steidel et al.||2010[|Rudie et al.||20 1 2) and 

provides important constraints for theoretical models which seek 


to determine the origin of gas in present-day galaxies l Keres et al. 


cosmological galaxy simulations remain systematically too high 
(e.g. |Keres et al.|2009| l. 

As resolving the relevant physical processes self-consistently 
has proven difficult, simulators have turned to physically-motivated 
subgrid recipes that allow the supernova energy from massive stars 
to couple effectively with the surrounding ISM and drive the galac¬ 
tic wind. In Smoothed Particle Hydrodynamics (SPH) simulations, 
the subgrid models which have been used can generally be divided 
into two classes. The first relies on injecting supernova energy as 
kinetic energy to selected SPH particles, and subsequently decou¬ 
pling the particles from hydrodynamic interactions to ensure their 


must also traverse the CGM, where it can run into and interact with 
outflows, before it can accrete onto galaxies. At later times, the 
CGM of massive halos contains a hot atmosphere, which can inter¬ 
fere with further accretion while also gradually cooling and supply¬ 
ing the galaxy. Recognizing the CGM as an environment where the 
IGM interfaces with the galaxy through this complex series of in¬ 
teractions has been a crucial step in properly framing the essential 
unsolved questions of galaxy formation jKeres et al.|2005| l. 

Deciphering the nature of galactic winds and the dynamic be¬ 
tween the ISM and CGM is particularly important in the context of 
regulating gas consumption during galaxy formation. Absent reg¬ 
ulation through strong stellar feedback, the known universal abun¬ 
dance of baryons in the ACDM framework would suggest that there 
should be many more galaxies with high stellar mass than what is 
observed (e.g. |White & Frenk|l99l| . Even when adding together 
the mass of all stars to the mass of the ISM, observed galaxies gen¬ 
erally do not appear to contain their share of available baryons in 
the cosmological context (liik ngila et al,| [l998; McGaugh et al,| 
|2010| >. Past efforts to account for the entire budget of "missing 
baryons" in the CGM and IGM have generally come up short, par¬ 
ticularly if it is required that the gas would be at least as hot as 
the virial temperature of L* galaxiesJBregmanJ2007 Anderson &| 


2005, Faucher-Giguere et al. 2010, 2011^ van de Voort et al. 2011, 

escape from the galaxy (e.g. Springel & Herne 

uist 20031 Oppen- 

Shen et al. 2013[ Ford et al. 2014; Faucher-Giguere et al. 2015). Fil- 

heimer & Dave|2006 

Oppenheimer et al. 2010; 

Vogelsberger et al. 

amentary infall from the high-redshift intergalactic medium (IGM) 

2013 see Dalla Vecc 

tia & Schaye 2008, Schaye et al. 2010 for 


|Bregman|2010|[Shull et al.|2012||Anderson et al.|2013 1 . However, 

the latest results from the HST COS Halos Survey hint that when 
taking into account ionization corrections, massive amounts of gas 
may be hidden in warm and cool clouds which together could ac¬ 
count for most, if not all of the missing baryons (Werk et al.|2014[ l, 
though the long-term stability of such a medium would likely re¬ 
quire it to be in a complex dynamical state. In galaxies less massive 
than L*, winds generated by stellar feedback are thought to sup¬ 
press star formation and eject gas from galaxies into the CGM and 
IGM jLarson|1974||White & Rees|1978||Dekel & Silk|1986| >. This 
hypothesis is used in all modern semi-analytical models of galaxy 
formation to match observed constraints on a wide range of galac¬ 
tic properties at various epochs (Benson et al. 2003; Bowe r~et al,| 
[2O06> . 

Hydrodynamic cosmological simulations promise to provide 
a direct, versatile, and comprehensive way of studying galaxy for¬ 
mation, but in practice, properly implementing realistic stellar feed¬ 
back physics that generates galactic winds has proven to be a sig¬ 
nificant challenge. Early efforts encountered the so-called "over- 
cooling” problem, where thermal energy injection from supernovae 
dissipated quickly (j Katz| | 1992[>, before it could generate a strong 
dynamical effect within the ISM (see Kim & Ostriker ( 20l5]> for 
a modern interpretation; see also |Martizzi et al,| ( |2015 ). Without 
strong dynamical effects from SNe that could generate galactic out¬ 
flows, the star formation rates and stellar masses found in many 


kinetic methods that do not decouple). The second approach uses 
the so-called "blastwave" model to account for the unresolved ex¬ 
pansion phase of the supernova remnant following the explosion 
(e.g. |Thacker & Couchman|2000||Governato et al.|2007[|Stinson 


et al. 1120061 IGuedes et al.||20 1 1 1 |Stinson et al.||20131 |Shen et~ 


2014[ Christensen et al.[|2015} . In practice, this means that after 
the thermal energy is injected into surrounding particles, cooling is 
disabled for a fixed length of time (other groups e.g.|Dalla Vecchia| 
|& Schaye|2012||Schaye et al.|2015| do not explicitly disable cool¬ 
ing but instead bottle up thermal energy until it is guaranteed to 
produce a blastwave). Simulations using either Adaptive Mesh Re¬ 
finement (fAgertz et al. 20 1 1 f or Moving Mesh (Vogelsberger et al. 
|2014) that have successfully run to z = 0 have had to rely on similar 
subgrid models. 

These subgrid approaches have, by construction, proven suc¬ 
cessful remedies for the overcooling problem, and have been used 
to produce galaxies with properties that are consistent with a vari¬ 
ety of observational constraints of late-type galaxies at z = 0. How¬ 
ever, these simulations have limited predictive power for proper¬ 
ties of the CGM. In large-volume cosmological simulations where 
star-forming regions are not well resolved (e.g. |Ford et ah]|2014[ 
|Vogelsberger et al.|2014} . the sub-resolution model encapsulates a 
significant portion of the galaxy. Galactic winds from these coarse 
(~ 1 kpc) regions are launched with unresolved phase structure, pre¬ 
determined uniform velocities, and pre-determined mass-loading 
efficiencies. Simulations that do resolve star-forming regions but 
employ the blastwave model (e.g. Stinson et al .|[2013| ) may over¬ 
heat the galactic wind material artificially as a consequence of tem¬ 
porarily disabled cooling (Agertz et al. 2013 Ma rtizzi et al.|2015) . 
Gas that has joined the CGM through galactic winds generated by 
these prescriptions may arrive with unrealistic properties. 

In order for simulations to play a role in improving our under¬ 
standing of the formation and dynamics of the CGM, particularly 
given the complex, multi-phase picture emerging from the latest 
observations ( [Tumlinson et al.|201 l[|Werk et al.|2014fr , the level of 
detail and sophistication in stellar feedback models must improve. 
In this work, we analyze the outflowing (and infalling) gas seen 
in the galaxies and CGM of the Feedback In Realistic Environ¬ 
ments (FIRE) simulation.^ first presented in Hopkins et al. (2014). 
Unlike the subgrid recipes which involve kinetically prescribed de¬ 
coupled winds and cooling-suppressed blastwaves, the FIRE sim¬ 
ulations solve the "overcooling" problem by explicitly modeling 
the radiation pressure, stellar winds and ionizing feedback from 
young stars as taken directly from the population synthesis code 
Starburst99 jLeitherer et al.|1999). These "early feedback" mecha- 
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nisms act before SNe, heating and stirring the surrounding interstel¬ 
lar medium which is necessary to match conditions in star-forming 
regions such as Ca rina (|H arper-Clark & M urray|2009} and 30 Dor 
(Lopez et al.|2011 ; Pellegrini et al. 2011} . SNe are implemented by 
taking into account their energy and momentum input. When the 
cooling radius of SNe is resolved, SN energy injected is free to ex¬ 
pand and generate momentum in the ISM before too much energy 
is radiated away. When this scale is poorly resolved, momentum ac¬ 
cumulation from SN remnant evolution below the resolution scale 
is added to the surrounding gas. This model is physically realistic 
when it is applied on the scale of giant molecular clouds, meaning 
that a resolution of several to tens of parsecs is required. The phys¬ 
ical feedback implementation in FIRE successfully regulates mass 
accumulation in galaxies and provides a physical explanation for 
the inefficiency of star formation in galactic disks ( jKennicutt[1983| 
|1998||Genzel et al.|2010| >. We stress that we allow hydrodynamical 
interactions and cooling of all gas at all times, unlike in typical sub¬ 
grid models. This is critical to make meaningful predictions for the 
phase structure of circumgalactic gas. 

This feedback model has been developed and tested in ideal¬ 
ized galaxy simulations l Hopkins et al.|20fT| |2012c|b|a[ see also 
|Agertz et al.|2013} . In |Hopkins et al. < 2014) (H14) the model has 
been extended and appropriately modified to be applied in a cosmo¬ 
logical framework. These simulations were able to reproduce the 
correct relationship between stellar mass and halo mass at a variety 
of epochs (see also |Agertz & Kravtsov|2015a|b| |Trujillo-Gomez| 
|et al.| |2015) and provide a way to explain the under-abundance 
of baryons in galaxies, relative to the cosmic mean predicted by 
ACDM (Fukugita et al.|[T998} |McGaugh et al.||2010} . Other suc¬ 
cesses of the simulations include reproducing HI covering fractions 
around Lyman-break galaxies (Faucher-Giguere et al.|2015} , mass- 
metallicity relations for both gaseous and stellar components of 
galaxies fMa et al.|2015} , and dark matter profiles that ease tension 
between ACDM and observed density structure of galaxies (Chan| 
|et al.|2015|[Qnorbe et al.|2tTi3) . 

Using these simulations, we explore the degree to which this 
stellar feedback model drives winds through the galaxy and the 
halo, their kinematics and phases, and the ways that these out¬ 
flows affect the long-term fate of the galaxy and the circumgalac¬ 
tic medium. We explore the relationship between gaseous infall, 
star formation, and outflow rates as measured in simulations. We 
will also attempt to bridge the gap between the methods typically 
employed by observers and simulators, and offer a framework for 
future comparisons between simulations employing various tech¬ 
niques to measure outflow rates. 

The present work will focus on a measurement of the mass¬ 
loading factor rj, which quantifies the ratio of the outflow rate to 
the star formation rate. We explore the nuances associated with 
measuring in simulations, and provide formulae for determining 
mass-loading factor and characteristic wind velocity as a function 
of halo properties and redshift. We also estimate the amount of gas 
each halo can deposit into its own CGM and how much of it leaves 
the virial radius and joins the IGM. The results of this analysis will 
be directly applicable to semi-analytic models and large-volume 
cosmological simulations. Phase structure, kinematics, and metal- 
licity evolution of winds will be explored in a companion paper. 

Section [2] will give a brief overview of the FIRE simulations. 
In Section[3] we discuss our method for measuring inflow and out¬ 
flow rates. We present analysis of inflow and outflow rates over 
time for a variety of halos in Section [4] We provide measurements 
of the mass-loading factor, r/, as well as wind velocities in Section 
[5] We discuss the long-term impacts of galactic outflows on the 


evolution of galaxies and the CGM in Section[6]and summarize our 
conclusions in Section[7| 


2 SIMULATIONS & ANALYSIS 

The simulations described in this work were performed with the 
pressure-entropy formulation of SPH as described in |Hopkins| 
(2013} , along with additional modifications detailed in H14. This 
formulation of SPH mitigates issues discovered in previous incar¬ 
nations around contact discontinuities and due to artificial viscosity. 
These fixes are particularly crucial for detailed study of the CGM, 
as the numerical artifacts of some older SPH codes produced dense 
clumps in the halo, and underestimated cooling rates from the hot 
halo atmosphere (Keres et al. 2012; Hayward et al. 2014; Keres et 
al. in prep). 

The simulations employ the zoom-in technique to follow se¬ 
lect Lagrangian regions around the halos of interest in consider¬ 
ably larger cosmological boxes. This technique limits our capac¬ 
ity to study the distant IGM, but captures the essential physics 
in circumgalactic gas flows, and gives us excellent resolution in 
the galaxy itself, including fully resolved structures in star form¬ 
ing regions. This enables us to use star forming threshold density 
of 10 — 100cm -3 for all of the galaxies analyzed here and imple¬ 
ment stellar feedback within giant molecular clouds. We analyze 
the same suite of simulations that was described in H14, many 
of which had initial conditions that were drawn from the AGORA 
suite (Kim et al.|2014} , which samples halos in decadal increments 
in mass from IO^Mq to 10 l ’Mq at z = 0. Particular emphasis 
is placed on halos that are in the L* range with mass ~10 i2 Mq 
at z=0 - the AGORA sample provides two such initial conditions 
(ml2q and ml2i), which we supplement with an additional set of 
L*-progenitor initial conditions (ml2v) studied extensively in prior 
work (Keres & Hernquist|2009| |Faucher-Giguere & Keres|2011} . 
The isolated dwarf galaxies m09 and mlO are extensively described 
in |Onorbe et al.|(2015} . To complement this sample, we also use the 
suite of simulations from |Faucher-Giguere et al.)(2015} , which only 
ran until z = 2, but provide a sample of relatively high-mass halos 
(up to 10 12 Mq at z = 2). We refer to this sample as the z2h sam¬ 
ple. Although the zoom-in technique generally focuses on a single 
“main” halo at z = 0, the multitude of progenitor halos at high red- 
shift were also analyzed, and included in the results presented in 
Section[5] We exclude ml3 (described in H14) due to concern that 
the resolution is too low to properly capture galactic wind genera¬ 
tion. 

The particle masses, softenings, and cosmological parameters 
employed in each run are give in Table [I] wh ich is in part repro¬ 
duced from H14, and |Faucher-Giguere et al.|(2015} . We also in¬ 
clude stellar and halo masses at z = 2 and z = 0 within R v u- for each 
halo. Note that a few values for halo mass differ slightly from the 
other papers due to different methods of halo finding. 

We use the public distribution of Amiga Halo Finder (AHF) 
(Knollmann & Knebe||2009} to determine halo centers and mea¬ 
sure central velocities using the adaptive mesh hierarchy method. 
We take the definition of virial overdensity from |Bryan & Norman| 
(l998j ), which evolves with redshift. We connect “main” halos to 
their progenitor halos using AHF’s merger tree module. We man¬ 
date that the physical value of the virial radius must stay constant or 
increase monotonically with time, to ensure relatively smooth vari¬ 
ation. We restrict most of our analysis to central galaxies, where 
this monotonic growth is expected, rather than satellite galaxies. 
For satellites, it is considerably more complicated to define quan- 
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Table 1. Simulation Initial Conditions 


Name 

M h {z = 0) 

M h (z = 2) 

M„(z = 0) 

M,(z = 2) 

m b 

£b 

m dm 

^dm 

Merger 


[%1 

[M r ] 

[M r ] 

[Mpl 

[M Q ] 

[pc] 

[%] 

[pc] 

History 

m09 

2.5e9 

1.3e9 

4.6e4 

4.1e4 

2.6e2 

1.4 

1.3e3 

30 

normal 

mlO 

7.8e9 

3.8e9 

2.3e6 

1.7e6 

2.6e2 

3.0 

1.3e3 

30 

normal 

mil 

1.4el 1 

3.8elO 

2.3e9 

3.4e8 

7.1e3 

7.0 

3.5e4 

70 

quiescent 

ml2v 

6.3el 1 

2.0el 1 

2.8el0 

2.3e9 

3.9e4 

10 

2.0e5 

140 

violent 

ml2q 

1.2el2 

5. lei 1 

2.2el0 

7.0e9 

7.1e3 

10 

2.8e5 

140 

late merger 

ml2i 

1.1 e 1 2 

2.7el 1 

6.1el0 

3.9e9 

5.0e4 

14 

2.8e5 

140 

normal 

z2h350 


7.9el 1 


9.0e9 

5.9e4 

9 

2.9e5 

143 

normal 

z2h400 

- 

7.9el 1 

- 

7.0e9 

5.9e4 

9 

2.9e5 

143 

quiescent 

z2h450 

- 

8.7el 1 

- 

1.3el0 

5.9e4 

9 

2.9e5 

143 

normal 

z2h506 

- 

1.2el2 

- 

1.8el() 

5.9e4 

9 

2.9e5 

143 

violent 

z2h550 

- 

1.9el 1 

- 

4.4e9 

5.9e4 

9 

2.9e5 

143 

quiescent 

z2h600 

- 

6.7el 1 

- 

1.7el0 

5.9e4 

9 

2.9e5 

143 

violent 

z2h650 

- 

4.0el 1 

- 

6.6e9 

5.9e4 

9 

2.9e5 

143 

normal 

z2h830 

- 

5.4el 1 

- 

1.4el0 

5.9e4 

9 

2.9e5 

143 

normal 


Parameters describing the initial conditions for our simulations (units are physical): 

(1) Name: Simulation designation. 

(2) My, (z = 0): Mass of the z = 0 '‘main" halo (most massive halo in the high-resolution region). 

(3) M/, (z = 2): Mass of the z = 2 "main" halo. 

(4) M* (z = 0): Stellar mass of the z = 0 "main” halo. 

(5) M* (z = 2): Stellar mass of the z = 2 "main” halo. 

(4) mt>: Initial baryonic (gas and star) particle mass in the high-resolution region, in our highest-resolution 
simulations. 

(5) ej: Minimum baryonic force softening (minimum SPH smoothing lengths are comparable or smaller). 

(6) m ( i , n : Dark matter particle mass in the high-resolution region, in our highest-resolution simulations. 

(7) z L h „: Minimum dark matter force softening (fixed in physical units at all redshifts). 


titles like the virial radius and the circular velocity, which we use 
to compare halos in Section[5] Furthermore, satellites are relatively 
poorly resolved in our simulations compared to zoomed-in isolated 
galaxies of the same mass. We have verified that satellites con¬ 
tribute relatively little star formation and CGM outflows compared 
to central galaxies. 

Star formation rates are computed by adding up the mass of 
all stellar particles that have been formed since the previous snap¬ 
shot, and then dividing by the time interval between the two snap¬ 
shots. To account for mass loss between the exact formation time of 
each stellar particle and the snapshot time, we boost the measured 
star formation rate by 15%, an approximation based on our typi¬ 
cal snapshot interval (-50 Myr) and the expected mass loss due to 
stellar evolution from a |Kroupa] ( |2001j l IMF. In this work, we con¬ 
sider only star formation within 0.2R v ir when comparing to outflow 
rates from the halo center. The star formation rate outside 0.2R v ir is 
also tracked in case it is relevant for anomalous outflow episodes. 
We compute star formation rates on a fine temporal resolution scale 
equivalent to the snapshot spacing. 

We quantify mass outflow rates from galaxies using the mass¬ 
loading factor, rj: 


where M ou t is the gas outflow rate and M* is the star formation rate. 
This quantity has been employed in both observational studies and 
in simulations to describe the strength of galactic winds. 

In |Hopkins et al.|(2012b) , the mass-loading factors were mea¬ 
sured in four isolated galaxy simulations, which featured a prelim¬ 
inary version of the feedback prescription employed in the FIRE 


simulations. In the present work, we aim for a more thorough in¬ 
vestigation of outflows in the CGM in a cosmological framework. 
Using the large suite of simulations, our data provide enough ha¬ 
los to produce statistically meaningful fits for various properties of 
outflows for halos at various epochs throughout cosmic time, which 
can be applied directly in semi-analytic models and larger simula¬ 
tions which cannot resolve the sources of feedback directly. We fo¬ 
cus on measuring a net flux at various points in the halo, including 
a measurement of how much gas leaves the virial radius. 

We also explore the nuances in measuring mass-loading fac¬ 
tors using a variety of methods, and show how the derived values 
are sensitive to spatial and temporal factors. The next section will 
describe the methodologies we have adopted for computing out¬ 
flow rates. Section 0 will demonstrate how these outflow rates are 
correlated with star formation in a way that can be quantified with 
a mass-loading factor r). Section[5]provides measurements of r], as 
well as wind velocity, for our large sample of halos. We also pro¬ 
vide fits that can be used in semi-analytical models or large-scale 
simulations, as well as to guide observational studies. 


3 MEASURING OUTFLOW RATES 

Using the Lagrangian property of SPH particles, it is possible to 
track parcels of gas directly as they drift through the simulated vol¬ 
ume (e.g. |Keres et al.|2005[ |Oppenheimer et al.|2010||Ford et ah| 
|2013) . Though the SPH particles employed in our simulations cer¬ 
tainly can be analyzed with particle tracking, we elect to use kine¬ 
matics to distinguish outflowing and infalling gas. We find this to 
be sufficient for determining instantaneous flow rates. Whether a 
particle is considered an outflow or an inflow is simply determined 
by the velocity in the radial direction relative to the halo center. 
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In other words, an SPH particle is counted as an outflow if at the 
epoch of the snapshot, it has a velocity v relative to the halo center 
such that: 


where r is the position vector from the halo center and |r| is the dis¬ 
tance from the halo center, and v cuf is the threshold velocity. Using 
specific values for v cut could allow a direct comparison with obser¬ 
vations at a given resolution in velocity, or where a known velocity 
dispersion in gas unrelated to the galactic wind can contribute to 
the rates. We use the same threshold to identify infalling gas, only 
reversing the sign of the velocity vector. Infall rates are given with 
a negative sign. 

We have experimented with a variety of values for this thresh¬ 
old, but we will primarily present results that use a simple v cu , = 0 
threshold to ensure that we account for a complete census of out¬ 
flowing and infalling gas. This simplifies comparison to and imple¬ 
mentation of our formulas in semi-analytical and subgrid models. 
However, in the interest of a more meaningful observational com¬ 
parison, and to quantify the amount of unbound gas, we also pro¬ 
vide measurements using v cuf = aw, the 1-D velocity dispersion of 
the halo. More information is found in Appendices |A| 

Once outflowing particles are identified, we calculate outflow 
rates in Mq yr -1 by computing the instantaneous mass flux through 
a thin spherical shell, in a manner similar to |Faucher-Giguere et al.| 
( [201 lj ). For simplicity, we divide the halo into finite spherical shells 
of thickness dL. The outflow rate is then computed with the follow¬ 
ing formula: 

= y> rr m sPH/dL. (3) 

at \r\ 

This value is summed up for all particles within the shell. 
When discussing flux measurements with this approach, we denote 
shells by giving their median halocentric distance. For example, 
the shell that encapsulates 0.2-0.3^ v ir is called the 0.257? v ir shell. 
For the analysis presented below, we use dL = 0. lfGir- We treat the 
innermost shell, 0 < r < 0. li? v ir, as the shell that contains the dense 
galactic ISM, while all outer shells are considered CGM halo gas. 
This analysis relies on the assumption that the outflows are gener¬ 
ated in this inner region of the galaxy by star formation that is also 
confined to the center. Alternate methods that explicitly treat ge¬ 
ometry and filling factors of outflows and inflows will be explored 
in future work. 

Our ability to do kinematic analysis relies on the time-spacing 
between snapshots being sufficiently short to capture each individ¬ 
ual star-forming and outflow event. Our typical interval between 
snapshots is between 10-40 Myr at z > 1 and 20-70 Myr at z < 1, 
which is sufficient for this analysis. In Appendix[A]we demonstrate 
that the outflow rates derived through this Instantaneous Mass Flux 
method are comparable to those computed using a comparable par¬ 
ticle tracking method that explicitly accounts for material crossing 
the shell, and determines the flux simply by accounting for all of 
this material and dividing by the time interval (Interface Crossing 
method). 


4 TRACING OUTFLOWS OVER TIME 

Our method for measuring outflow rates requires that outflowing 
material has a detectable kinematic signature in the snapshots of 


our simulation. This method would not work if the interval between 
snapshots were long enough such that gas particles launched as 
winds during the interval would slow down and randomize their ve¬ 
locity vectors prior to the next output. To demonstrate that we avoid 
this pitfall, we must show that particles marked as outflows are 
both sufficiently kinematically distinct from "static" gas at the same 
galactocentric distance, and that "outflowing material" as marked 
by the detection method really does move outward, while inflow¬ 
ing material also has the correct kinematic and spatial structure. 

In Figure [T] we show gas as binned by radial velocity (v ra d) 
and weighted by the "flow rate" from Equation [3] at four epochs. 
Alongside each of these flow rate distribution diagrams, we show 
the surface density of outflowing and infalling gas, in a 2-D projec¬ 
tion of the halo. The halo shown is one of the two main progeni¬ 
tors of nil2v immediately following a major burst of star formation 
(peaking at ~12Mq yr -1 ) atz « 3. The first and second rows of Fig- 
ure[l]show the halo at an epoch in the early stage of a major outflow 
episode, when outflowing material is still confined to the ISM (first 
row), and approximately 50 Myr later (second row), when material 
has reached the CGM. In the early stage of the burst, the only mass 
flow seen in the CGM is a small gas inflow with v ra d ~ aw- The fil¬ 
ament that is feeding the galaxy with fresh material remains intact 
at this stage. This changes drastically in the second panel, as a vast 
swath of outflowing material has reached the inner CGM. This can 
be seen in the spatial 2D projections, as the 0.257? v ir shell now con¬ 
tains several clumps of high-column density outflowing material, 
while the filament has been disrupted in these inner regions. The 
Vmd histogram shows that a large amount of material now occupies 
a broad distribution centered at v m d ~ 250km/s, significantly offset 
from a id- The fastest material has already streamed past that shell, 
as can be faintly seen in the v m d flow-rate diagrams in the 0.57? v ir 
shell. 

In the third and fourth rows of Figure |Tj we see the continu¬ 
ation of the outflow episode at -100 Myr and -150 Myr after the 
outflow began. In the v ra d flow rate distribution, gas appears to be 
streaming into the 0.57? V ir shell in the third row, and the 1.07? v ir shell 
in the bottom row. The velocity flow rates have the same charac¬ 
teristic broad distribution as was seen in the 0.257Air shell earlier, 
though with lower outflow rates, suggesting that a large fraction 
of the material does not escape the virial radius of the galaxy and 
either remains in the CGM or recycles. Infall rates over the inter¬ 
val of time shown remain low, but there is indication that a net 
flow of material is starting to trickle back into the 0.25/? v ir shell, 
while the filamentary inflow geometry is restored at the last epoch 
considered. The characteristic velocity distribution of inflow never 
reaches the high speed of outflows. 

To best understand the nature and prevalence of outflows 
in the CGM, and to connect them with galactic evolution, 
measurements over a time interval much longer than a single 
outflow episode is required. For each snapshot of each halo in 
our simulation suite, we calculate the star formation rate, outflow 
rate and inflow rate at various radii and plot them as a function 
of cosmic time. Inflow and outflow rates are measured using the 
Instantaneous Mass Flux method (Equation [3} with v cu , = Okm/s 
unless otherwise stated. The first of such time evolution plots 
is shown in Figure [2] This Figure follows the progenitor of one 
of the L* galaxies, ml2v, during the early stages of formation, 
4.5 > z > 2.0. The extremely bursty nature of star formation and 
"gusty," as opposed to steady, nature of galactic outflow rates are 
the most obvious results, and we shape the remainder of analysis in 
this paper based on this behavior. We enumerate several interesting 
features of the plot: 
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Figure 1. Left: histogram showing the distribution of halo gas radial velocity, weighted by the flow rate before (first row) and during (second, third, and fourth 
rows) a major galactic wind flows through the CGM at z ~ 3 in an L*-progenitor halo (ml2v). Vertical dotted lines represent v ra d = d=cr 1Z) for the halo. 
Histograms are shown in three thin shells of width dL = 0.1/? v ir> centered at 0.25R v ir (black, solid), 0.57? v ir (blue, dashed), and 1.0/? v ir (red, dash-dot). Each 
row represents a snapshot at a particular epoch, with a spacing of ~50 Myr between epochs. Center: corresponding to the same epoch as each radial velocity 
histogram, the mass surface density per 2-D spatial bin computed on a 300 x 300 grid between —R w j r < x,y < R wlT . Only particles with \v ra d\ > &\d are 
shown. Right: mass surface density of infalling gas computed with the same methodology as the outflowing particles in the central panel. White dotted circles 
denote 0.25R v ir and 0.5R v ir- The thick white circle denotes R v ir- 
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Figure 2. Star formation, outflow, and inflow vs. cosmic time in one of 
the L* progenitor halos, ml2v, for 4.5 > z > 2.0. Star formation is shown 
in black, outflows rates in red at 0.25/?™ (dashed) and 1.02? v ; r (dotted). 
Infall rates are shown for the same two halo-centric radii in blue, and are 
represented as negative flux. The star formation rate has been rescaled by a 
factor of 10 to be visible on the same linear scale. Outflow rates in the CGM 
are variable, and peak after bursts of star formation, though with a temporal 
offset. Inflow rates are much more steady, but can be disrupted temporarily 
by outflows. 


1. The bulk of star formation occurs in bursts. The galaxy 
spends a significant fraction of time with near-zero star formation 
rates. 

2. Bursts of star formation are always followed by outflows at 
0.25/?vir, and the outflows rates are typically higher than the star 
formation rates. 

3. Bursts of star formation are preceded by inflow at 0.25/?™. 

4. Outflows through 0.25/? v ir are temporally offset from star for¬ 
mation. They would be observed as coincident within the timespan 
(100 Myr) of some star formation diagnostics (e.g. UV excess), 
but not within the timespan of others (e.g. Ha). We find outflow 
rates in the inner CGM at high redshift to be most correlated to 
star formation at a time dr fts 60Myr after each burst, although the 
scatter about this timescale is large (see Appendix |B 4 . 

5. Outflows through 1.0/?™ are even more temporally offset from 
star formation, sometimes by a longer period (-250 Myr) than 
common star formation diagnostics pCennicutt & Evans|201 2\ . 

6. Infall in the inner halo (0.25/? v ir) is temporarily suppressed when 
outflow rates are high. This suggests that inflows can geometrically 


overlap with outflows and interact with them (see also Oppen- 


heimer et al ,|20 1 (J) |van de Voort et al.||20l 1| Faucher-Giguere & 


Keres [2011} Faucher-Giguere et al.|2011 |2015}. 

7. Outflows at 0.25/?™ do not always reach 1.0/?™. Even when 
the outflows do reach 1.0/?™, the outflow rates are typically not as 
high, and instead persist in the shell for a longer period of time, 
suggesting a slower flow rate in these outer regions. 

8. Outflow, infall, and star formation rates are all highly variable. 
Throughout the 4.5 > z > 2.0 interval shown, this galaxy never 

switches into a steady mode of star formation. We will later show 
that more massive galaxies may have continuous star formation, 
particularly at lower redshift. 


Taken together, these features suggest a highly time-variable 
physical state of halo gas. Even the inner region of the halo (the 
galaxy) is prone to instability. Outflows are typically launched from 
within the galaxy, pushing out and dispersing the ISM. We find 
that during the 4.0 > z> 2.0 interval represented in Figure [2] the 
amount of dense gas within /?™ with n > lcm ~ 3 only exceeds 
1 O 7 8 M 0 for short periods of time, corresponding to the snapshots 
when the star formation rate is high, and the snapshots immediately 
prior. In fact, the galaxy can be completely devoid of this dense ISM 
for several hundred Myrs following a major outflow episode. 

The inflow rates at 0.25/? V ir do not vary as drastically as out¬ 
flow rates. The inflow rate at 1.0/?™ is even more steady, only 
changing significantly as the system approaches a merger event 
at z ~ 2. The integrated mass flux going through 0.25/? V i r as in¬ 
flow slightly exceeds the integrated mass flux going out. This trend 
is even more pronounced when considering the integrated flux at 
1.0/? v ir, which means that the overall amount of baryons contained 
within /?™ is increasing. The implications of the outflows and in¬ 
fall for the overall baryon content of galaxies will be explored in 
Section f6.2l 

We stress that this galaxy has a halo mass of approximately 
2 x 10 1 *Mq and a stellar mass of 2.4 x 10 9 Mq at z=2, which 
is somewhat less than the mean Lyman-break galaxy halo mass 
(-IO^Mq) and therefore would fall below the selection threshold 
of many surveys ( jSteidel et al.||2003) . Therefore, the time vari¬ 
ability seen in this object is not necessarily representative of cur¬ 
rent LBG samples. H14 show the star formation history of the 
same galaxy smoothed over 100 Myr timescale, which more closely 
corresponds to commonly used observational star formation rate 
indicators. They find that the measured star formation rates are 
broadly consistent with the star-forming main sequence jNoeske] 
[et al.|2007| l. 

4.1 Time evolution of outflows 

To follow the evolution of SFR, outflows, and inflows, we divide 
the redshift range being considered in our analysis, 4.0 > z > 0.0, 
into three increments. The first, "high-z" interval is 4.0 > z > 2.0, 
where our dataset is the most complete due to the inclusion of the 
z2h sample and large number of progenitors and satellite progeni¬ 
tors in the zoom-in region. The next division is somewhat arbitrary, 
but is justified based on the general nature of outflows and star for¬ 
mation in the L*-progenitor halos, which we will explore in subse¬ 
quent sections. We use 2.0 > z > 0.5 as the intermediate redshift, 
or "med-z" interval and 0.5 > z > 0 as "low-z" interval. We choose 
not to consider epochs at z > 4.0 to maintain a high standard for 
resolution elements per galaxy, and to confine analysis to typically 
observable redshifts. 

During the high-z interval, all halos in our sample undergo 
bursty star formation and gusty outflows, interspersed with peri¬ 
ods of quiescence as seen in Figure [2] We plot several of the halos 
spanning a wide range of masses in Figure[3] Though the scales of 
activity are different, it is clear that high time variability is ubiqui¬ 
tous. This is generally consistent with observations of highly active 
low-mass star-forming galaxies at high redshift (e.g.|Maseda et al.] 
|20l4) . 

The relative efficiency of galactic wind generation depends on 
the mass of the halo. At one end of the spectrum, dwarf galaxies 
like mlO (upper left panel of Figure [3j and m09 (not shown) need 
only form a small number of stars to expel the entire ISM and tem¬ 
porarily shut off star formation altogether. At the other end of the 
spectrum, we consider the most massive galaxy in the z2h sample, 
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Figure 3. Same as Figure[2] but for four other halos (upper-left: mlO, lower-left: mil, upper-right: ml2i, lower-right: z2h506). In all cases, the star formation 
rate is rescaled by a factor of 10 to appear visible on the linear scale. While the scales needed to show SFR and outflow rates for each figure differ by orders of 
magnitude, the pattern of bursty star formation and gusty outflows appears is generic. In the most massive halos, star formation does not shut off after z = 3. 


z2h506 (lower right of Figure [ 3 ] which has a mass of 1O I2 M0 at 
z = 2). In this halo, star formation, inflow, and outflow seem to co¬ 
exist in a continuous, albeit highly variable state. The rest of the 
halos generally behave in a manner between these two extremes, 
and that behavior can largely be predicted by halo mass. This is 
expected, as halo mass is a good proxy for the depth of the grav¬ 
itational potential well at the center of the galaxy, and for the rate 
of gas accretion from the IGM, both of which significantly impede 
galactic outflow launching and propagation. It is interesting that in 
the most massive halos in our sample at z = 2 at around 10 i2 Mq, 
shutdown of star formation after bursts is not complete. While the 
majority of stars form in bursts, some level of continuous star for¬ 
mation rate is also being established. 

In the med-z regime, which is shown in the top panel of Figure 
[4] the L*-progenitor ml2i continues the pattern of bursty star for¬ 
mation and gusty outflows seen in the high-z regime, especially at 
z > 1. However, the outflow rate gradually decreases relative to the 
star formation rate suggesting a decreasing t;. Furthermore star for¬ 
mation enters a continuous, non-bursty mode starting at z « 0.75. 
This halo behaves significantly differently than the previously men¬ 


tioned massive high-redshift z2h506 halo, though the mass of ml2i 
at z = 0.5 is comparable to the mass of z2h506 at z = 2. While 
z2h506 drove powerful, gusty outflows and simultaneously main¬ 
tained ceaseless star formation at high redshift, ml2i has entered 
into a continuous star formation mode with weak galactic winds. 

Other L*-progenitors, ml2q, ml2v, and the massive dwarf 
mil also have reduced outflow rates relative to their star formation 
rates, but continue the bursty pattern of star formation and power¬ 
ful outflows through the CGM during the interval. One interesting 
feature of this regime which will be elaborated on in Section [6~2] 
is that despite a lower efficiency in wind generation, all of these 
halos appear to lose more material through outflows than they gain 
through inflows at /? v i r . This is due to lower inflow rates from the 
intergalactic medium, combined with a few isolated incidents of in¬ 
tense star formation which generates the largest increases of stellar 
mass in the history of the halos, followed by the most intense galac¬ 
tic superwinds ever seen in the halos. The dwarf galaxy mlO also 
continues to form stars and generate outflows and bursts, though 
they are relatively diminished compared to higher-redshift activity. 
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Figure 4. Same format as Figure [2] but showing the 2.0 > z> 0.5 and 
0.5 > z > 0 regime for the same halo (ml2i). In the low-z regime, we also 
show outflow rates measured with v cu1 = aw- Compare this figure to the 
upper-right panel of Figure [5] to see the full evolution of ml2i with red- 
shift. Outflows in the CGM become increasingly uncommon at low redshift, 
while star formation switches into a continuous mode. Note that unlike pre¬ 
vious figures, star formation rate is not rescaled by a factor of 10. In lower 
mass galaxies such as mlO and mil, the bursty nature of star formation and 
gusty outflows persist at low redshift. 


Unlike the other dwarf galaxies, m09 does not form a significant 
mass of stars near the galactic center after z = 2. 

In the low-z regime of ml2i, which is shown in the bottom 
panel of Figure]?] a qualitative shift is seen in the nature of star for¬ 
mation and outflows. Stars form at a moderate, steady rate, while 
outflow rates in the CGM are low and no longer gusty. Though 
they are not shown, the behaviors of ml2q and ml2v are both sim¬ 
ilar to that of ml2i at this epoch, but with lower star formation 
rates. The inflow rates to all halos gradually taper off during this 
epoch, which coincides to the thinning nature of cosmic filaments 
and longer cooling times of the halo gas, which slows down gas 


infall to the 

inner CGM and the galaxy (Keres et al. 2005 ; Dekel 

& Birnboim 

2006; Faucher-Giguere et al.J^Oll, van de Voort & 

Schaye 2012 

Nelson et al. 2013). We will explore the nature of the 


remaining infall in future work. 


We consider one explanation for the drop-off in gusty outflows 


at 0.25fi vi r below z < 0.5 for all three L*-progenitor halos: that gas 
is still being ejected from the disk, but no longer able to traverse the 
CGM. In massive (10 12 Mq) halos, 0.25R vlr represents a great phys¬ 
ical distance from the central galactic core, particularly at low red¬ 
shift (-50 kpc). To test this scenario we measure the flux at 0.05R v i r 
and find that while the outflow rates in the inner regions are greater, 
there is no indication that the outflow at either 0.05f? V ir or 0.25/? v ir 
represents a kinematically distinct component of gas. This is fur¬ 
ther evidenced by the drastically reduced values of outflow rates 
when they were measured with v cu t = <j\d, where aw ~ lOOkm/s 
for this halo at z = 0. The outflow rates using this threshold for 
both 0.05f? v ir or 0.25i? v ir are shown in Figure [4] The strongly di¬ 
minished flux measured with v cu , = aw indicates that the apparent 
“outflows” at late times can be explained by small random veloci¬ 
ties of the halo gas as well as the motion of the material associated 
with and stripped from the satellite galaxies moving through the 
halo. Although the flux is somewhat higher at 0.05i? v ir, this can be 
entirely explained by the late-time accumulation and random mo¬ 
tions of the ISM. In both cases apparent outflows are not connected 
to the star formation events. We have explicitly checked that both 
of these scenarios contribute to the measured outflow rates and will 
further discuss this point in Section [5~2] 

The drop-off in CGM outflows is likely caused by the deepen¬ 
ing gravitational potential well in the center of the halo, the sparse¬ 
ness of gas-rich mergers (that trigger bursts of star formation), less 
geometrically concentrated star formation, and the overall more 
continuous, less bursty mode of star formation. We briefly explore 
the changing behavior of star formation by considering the pres¬ 
ence and quantity of dense gas. In the high-redshift regime, the 
mass of ISM with number density n > 1cm -3 fluctuates rapidly - 
ISM mass builds up, precipitates a burst of star formation, and be¬ 
comes largely depleted after each burst. Specifically, the ratio of 
the mass of gas with number density n > lcm~ 3 to galaxy stellar 
mass reaches 0.2-1 at the epochs prior to and during each burst, 
but can fall below 0.01 afterwards, and remains low during quies¬ 
cent periods. At low redshift, as ml2i has significantly increased 
in mass and central density jChan et al.|20l5> , dense gas can per¬ 
sist. The ratio of dense ISM to stellar mass remains 0.1-0.2 at all 
times for z < 0.5. Star formation is clearly no longer regulated by 
the same gas dispersal and ejection process that is so effective at 
high redshift. 

Instead, ml2i appears to regulate star formation via a stable 
gaseous disk at low redshift, which can now persist for significant 
timespans, as merging activity has largely ceased for L* halos. The 
emergence of such stable gaseous disks at late times is consistent 
with observational studies jKassin et al.|2012) . Feedback from star 
formation contributes to the stability of this disk by stirring tur- 
bulence ^Thompson et al.|2005[ Ostriker & Shetty|20fT||Faucher-| 


|Giguere et al.|2013[|Hopkins et al,|2013[ l. Gas affected by stellar 

feedback at low redshift can be launched as a local fountain, but 
generally does not propagate into the CGM. Measuring the preva¬ 
lence of these local fountains and quantifying their effects on the 
galaxy requires sophisticated methodology and is left for future 
studies. In this stable disk, star formation no longer happens in a 
dense clump near the center, but is more extensively distributed. 
Extended gaseous disks are also seen in ml2v and ml2q, although 
they are considerably less massive, and the bulk of star formation 
is primarily confined to nuclear regions. 

We note that the dwarf galaxies in mil and mlO maintain a 
bursty star formation history and outflows at low redshift, though 
the bursts acquire a new characteristic feature. Instead of a single 
burst of star formation followed by a corresponding outflow in the 
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Table 2. Mass-loading factors 


Name 

Flux T 

Cross T 

Flux T (T\D 

FluxI 

Flux E 

m09 

560 

790 

240 

180 

350 

mlO 

120 

140 

70 

170 

210 

ml 1 

14 

21 

10 

17 

14 

ml2v 

9.3 

9.3 

6.8 

9.1 

11 

ml2i 

8.0 

10 

5.9 

11 

9.2 

ml2q 

8.3 

8.5 

6.7 

6.6 

6.8 

z2h350 

8.6 

8.8 

6.7 

7.9 

7.6 

z2h400 

6.6 

5.7 

5.1 

9.1 

6.2 

z2h450 

6.4 

5.9 

5.1 

7.4 

6.6 

z2h506 

6.5 

6.9 

4.5 

6.4 

4.9 

z2h550 

7.6 

7.7 

5.5 

11 

7.4 

z2h600 

6.0 

6.3 

4.6 

7.2 

6.7 

z2h650 

6.8 

8.8 

5.1 

7.1 

7.9 

z2h830 

4.8 

4.2 

3.6 

5.3 

3.6 


Method used to calculate rj at 4 > z > 2: 

(1) Flux T: Outflows measured by Instantaneous Mass Flux 
method for all gas with v m j > 0. // is computed by the 
integrating total star's formed and total outflow rate over the 
whole interval. These values of rj were used in Section [5~T| See 
Section[3]and Appendix [aTTT] 

(2) Cross T: Same as (1) but using Interface Crossing Method 
to calculate outflow rates over the interval. See Appendix |Al ,2| 
(3 Flux T cf\d : outflows rates over the interval are measured 
using the Instantaneous Mass Flux method, but only particles 
with v ra i > (J i/> are counted as outflows. See Appendix |A2| 

(4) Flux I: Using Instantaneous Mass Flux with v cu t = 0 to 
measure outflow rates over the high-z interval, rj is computed 
by finding correlation of star formation rate at each snapshot to 
outflow rate at a time dr later (the instant correlation method, 
see Appcridix |B 1 1 . 

(5) Flux E: Using Instantaneous Mass Flux with v cut = 0 to 
measure outflow rates over the high-z interval, rj is computed 
by integrating total stars formed and mass expelled over 
individual episodes, and then deriving the correlation 
parameters (the episodic correlation method, see Appendix|B2l. 


CGM, it now takes a series of small bursts interspersed with pe¬ 
riods of quiescence to drive a sufficiently powerful outflow. This 
could indicate that gas is ejected from the ISM, but with insuffi¬ 
cient velocity to kick it out to large radii, until finally there have 
been enough episodes that a threshold for velocity or temperature 
is reached, and the outflowing gas escapes to the CGM. In both 
cases, the magnitude of star formation and outflow activity is con¬ 
siderably weaker than at earlier epochs. 


5 MEASUREMENTS OF MASS LOADING IN GALACTIC 
WINDS 

We now provide measurements of the mass-loading factor, rj, for all 
halos in our simulations. We note that there are at least two distinct 
approaches which can be used to measure tj, each of which is ap¬ 
propriate for different purposes. The first approach is to make an in¬ 
stantaneous measurement of rj - a direct measurement of how much 
material is ejected for a particular star formation episode. This ap¬ 
proach is useful for observational comparisons to data where rj 
is inferred from the current state of outflowing gas. On the other 
hand, this measurement is subject to stochastic effects, as each out¬ 
flow episode is influenced by highly variable circumstances of the 
galaxy’s dynamical state. As such, there is significant systematic 


scatter associated with these measurement. Nonetheless, we have 
made measurements of ij using several variations of this approach, 
and provide detailed methodology and results in Appendix |B1| and 
Table [2] 

The second approach, which we will focus on here, is to mea¬ 
sure an average cumulative rj over a time interval for each individ¬ 
ual halo. This measurement provides a useful diagnostic that can be 
integrated to recover the gas mass ejected into the CGM and out of 
the halo for some portion of its evolution. The main caveat to this 
approach is that in order to get a statistically significant represen¬ 
tative sample of star formation and outflow events for each halo, a 
relatively broad redshift interval must be used. 

We focus here on average cumulative values of rj over vari¬ 
ous intervals for halos in our simulations. We choose three such 
intervals to divide the 4.0 > z > 0 evolutionary history, which were 
already introduced in Section |4~T| as high-z, med-z, and low-z. Each 
of these intervals typically provides at least -5 distinct episodes of 
outflows from which to make a measurement; however, the halo 
can gain a very significant amount of mass between the start of the 
interval and the end. For example, the mass of the "main” halo of 
ml2i grows by a factor of 10 between z = 4.0, the start of the high 
redshift interval, and z = 2, the endpoint. This growth affects the 
physical conditions that affect the launching of the galactic wind. 
For this reason, we give fits for 77 as a function of various physi- 
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Figure 5. Average mass-loading factor at 4.0 > z > 2.0 (top, black), 2.0 > z > 0.5 (middle, blue), and 0.5 > z > 0 (bottom, red) vs circular velocity (left) and 
halo mass (right) as they are at the midpoint of the redshift interval (zmed = 3, 1.25 and 0.25 for high, medium and low-z). Right side up triangles show the 
halos in the zoom in regions of m09, mlO, and mil. Upside down triangles, show ml2 halos, except at low-z, where their “main” halos are shown as upper 
limits (see text). Squares show z2h halos. Open symbols indicate halos that survived at least as long as the midpoint of the interval, but not until the end of the 
interval. For 77 as a function of v c , we provide a broken power law fit (dotted line) including a redshift evolution term (Equations 4 and 5), and evaluate it at 
z = 3 (top), z= 1.25 (middle), and z = 0.25 (bottom). Fits were generated with high and intermediate redshift data and extrapolated to low redshift. The 77 vs 
Mf, fit is derived from the 77 vs v c fit (Equations [b] and ^}. 
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Figure 6. Average mass-loading factor (vj) vs stellar mass {M t ), using the 
same symbol and color conventions as Figure [5] A single power law fit 
with no redshift dependence (Equation [8] dotted black line) describes the 
data well, except for massive halos at low redshifts, where outflows are 
diminished (red upper limits). 


cal quantities as measured at the midpoint of each interval in red- 
shift (Zmed = 3 for high-z, Zmed = 1-25 for med-z, and Zmed = 0.25 
for low-z). We have considered other choices for the representative 
redshift, such as the epoch when the cumulative time-integrated 
flux of ejected material in each halo reaches 50% of its final value, 
but found that our results were largely unchanged. Within each red¬ 
shift interval, we elect to use a single epoch for all halos to simplify 
interpretation. 

In the figures and fits provided in the sections below, we 
present r] as measured by the ratio of integrated outflow and 
star formation rates over the entire considered interval. Outflow 
rates themselves were measured with the Instantaneous Mass Flux 
method, and a radial velocity cut of v cut = 0 is used to define out¬ 
flows. We also provide Table [2] which shows average values of 7 ) 
for the “main” halos in each simulation at various epochs using var¬ 
ious measurement methods. All outflow rates were measured in the 
0.25i?vir shell. Section[53]shows how these measurements differ at 
various halo-centric radii. An alternative approach would be to in¬ 
stead use a shell at a fixed physical radius at all times (i.e. a few 
tens of kpc). However, using such a threshold would probe rather 
different spatial regions when applied to our dwarf galaxies (poten¬ 
tially outside R v ■„), and to our most massive halos (close to galactic 
edge). For now, we stick to using shells at a fixed fraction of /? V ir, 
as they can consistently be adapted to all halos at all epochs. 


5.1 Fits of ?/ for individual halos 

We start by considering the relationship between r/ and the halo 
circular velocity (iy = \JGMi,/R v i r ), which evolves more slowly 
with redshift than other halo properties (as previously mentioned, 
the halo mass of ml2i increases by a factor of -10 between z = 4 
and z = 2, while v c only increases by a factor of -2). We show the 
average value r) vs ty in the left panel of Figure [5] We can imme¬ 
diately see that halos with low v c and halos with high v c may be 
best described by different slopes. Our method for constructing the 


fit for 7 ) vs. v c is as follows: We divide the sample into two distinct 
populations, v c < 60km/s and v c > 60km/s. The choice to use 60 
km/s was determined to produce the lowest \ 2 statistics compared 
to other dividing points. We then "anchor" both fits to the approx¬ 
imate mean value of r/ at ty = 60km/s, which is computed by tak¬ 
ing the values of eta for all halos between 50km/s < rj < 70km/s 
and averaging them in log space. We use this anchor to fit a bro¬ 
ken power law for low v c and high v c halos. We find that the re¬ 
duced \ 2 statistic for the broken power law fit is significantly lower 
than the reduced x 2 statistic for a single power law fit. All fitting is 
done in log-log space, and the errors quoted come from the diago¬ 
nalized covariance matrix from the minimized residuals. First, we 
consider only our high-z sample, which is the most complete (black 
points in Figure[5]i. Taking measurements of v c at z = 3, halos with 

v c < 60km/s are fit by the relation 77 = 17 . For halos 

1.0 


' ( —tc _/ 

y 60km / s J 


with v c > 60km/s, we find 

Because the z2h sample ends at z = 2 and halos begin to co¬ 
alesce and grow to encompass the zoom-in regions, the number of 
halos that we are able to measure is smaller at z < 2. Although our 
data at lower redshift is relatively limited, we also provide a fit for 
the med-z (2.0 > z > 0.5) sample. We again take measurements of 
v c at the midpoint of the redshift interval, in this case z= 1.25, and 

-3.0 


use 60 km/s as a dividing point. We find r/ = 8 


^ 60km/s) 


for 


Vc < 60km/s, and r) = 8.8 ^ 60 ^ //s J for v c > 60km/s. Since the 
slopes measured for the high-z and med-z samples are similar in 
both Vc regimes, we combine the datasets and provide a single uni¬ 
fying fit including an additional term to account for redshift evolu¬ 
tion. For v c < 60km/s, we find: 


7) = 2.9(1 +z ) 1,3 


/ Vc 
\ 60km/s 


-3.2 


For halos with v c > 60km/s, we find 


(4) 


’- 2 - 9 , 1 +d ”(«»sv;) ■ ,5 > 

The derived errors are a = 0.2 and a = 0.07 for the low-v c 
and high-iy sides, respectively. The error on the redshift evolution 
is a = 0.1. The root-mean squared error around the fit is 0.1 dex. 
If we fit halos with v c < 60km/s and v c > 60km/s separately and 
allow the normalization to be another free parameter, the errors for 
the two slopes are larger by a factor of -1.5, and the normalizations 
of both fits have errors of 0.2 dex, suggesting a factor of -2 uncer¬ 
tainty. We stress that these fits apply only to the mass-loading factor 
as measured at 0.25/? v i r , and would likely be somewhat different if 
we instead measured rj at a fixed physical radius^Fits using rj com¬ 
puted only from outflows with v ra d > &id show similar slopes for 
both sides, but are normalized lower by -25%. 

We choose not to include the limited dataset we have from the 
low-z (z < 0.5) regime in the fitting procedure, but note the most 
significant outliers are two of the three “main” L*-progenitor halos. 
In fact, the outflows seen in the L*-progenitors at this epoch are 
probably not caused by star formation, so we therefore mark them 


2 We have also measured 77 at a fixed physical radius of 25 kpc for all 
halos at all redshifts, and verified that 7 / at fixed ty has redshift evolution 


consistent with 


ith Equations|4]and|3] 
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as upper limits on the figures. This discrepancy will be explored 
further in Section [5~2l 

One notable feature of Figure[5]is that the simulations plotted 
were run with different resolutions, yet low-mass halos from the 
z2h sample as well as those from the L*-progenitor runs (ml2v, 
ml2i, ml2q) are found to have similar values of 77 as the highly- 
resolved dwarf galaxies (m09, mlO, mil). This demonstrates a de¬ 
gree of convergence in our simulation sample. 

We convert the broken power law fit derived for v c into a fit 
for halo mass, Mi, by using the analytic relationship between the 
two. The result is shown in the right panel of Figure [5] This fit 
works just as well as the v c fit by construction, owing to the one-to- 
one correspondence between v c and Mi, as measured at a particular 
epoch. We only need to consider the redshift evolution of Mi, at the 
fixed value of v c = 60km/s, which we represent as Mi, 6 o. 


For Mi, < Mi, 60 , we find: 



77 = 2.9(1 +z ) 13 

(AAV L \ 

\M h 60 / 

( 6 ) 

For Mi, > Mi, 60 , it becomes: 



77 = 2.9(1 +z ) 13 ( 

' Ml, \ ~° 33 
\Mi,60 ) 

(7) 


It is straightforward to derive Mi,so as a function of redshift for 
any given cosmology and choice of definition of virial overdensity. 
We find Mi,m = 1.8 x 1O 1 O M 0 at z = 3, M h60 = 4.0 x 10 10 M© at 
Z = 1.25, and M h60 = 8.3 x 1O 1 O M 0 at z = 0.25. 

Figure [ 6 ] shows the relationship between rj and stellar mass, 
M,. Unlike the fitting method used for v c and Mi,, we use a sin¬ 
gle power law fit that describes 77 as a function of M*. We have 
confirmed that the reduced x 2 statistic for this fit is low, validat¬ 
ing our single power law approach. We again combine the data for 
the high-z and med-z samples, and include an additional term in 
the fitting function for a redshift dependence. The best-fit relation 
for the redshift dependence is r] oc (1 +z) 0 ' 02 , which is consistent 
with no dependence, given our errors. Hence, we present a redshift- 
independent fit for 77 as a function of M»: 


/ \ -0.35 

” = 3 ' 6 (lo%) ' <8) 

The error on the fit for the power-law indices is a = 0.02, 
suggesting a better fit than what was found for v c , though the scatter 
is still significant. The error on the normalization is 0.2 dex. The 
L*-progenitors at low redshift are again the most notable outliers 
from the fit, while the rest of the low-z data is somewhat better 
described than it was by the v c fit. 

We note that when a halo has multiple long-lived massive pro¬ 
genitors that coalesce by the end of the interval, we include each 
progenitor as a separate track in the figures and fits. We employ a 
few other criteria for inclusion of progenitor halos: 

1. ) They are either the “main” progenitor for the entire interval, or 
are detected as central (non-satellite) halos until the midpoint of the 
interval in redshift space (z = 3 for high-z z = 1.25 for med-z, and 
z = 0.25 for low-z.) 

2. ) Consist of at least 98% high-resolution dark matter (low con¬ 
tamination from low-resolution particles). We make an exception in 
the case of mil and ml2q, as these two runs feature low-resolution 
DM particles that are about as massive as hi-res DM particles from 
other runs (ml2v and ml2i). Each consists of more than 98% high- 
resolution DM particles within 0.1 R v i r . 

3. ) Contain at least 50,000 high-resolution dark matter particles at 
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the end of the interval, or at the last epoch at which they are counted 
as isolated halos. 

4.) Form at least 50 new stellar particles over the interval. 

Our statistics at the high-z interval are sufficient to obtain a 
statistically meaningful fit even if we consider only halos that sur¬ 
vive as their own “main” progenitors to z = 2 , but inclusion of the 
other progenitors does not appear to alter the fit. We mark halos 
that did not survive until the end of the given interval as open sym¬ 
bols on the figures. Although simulations of different resolutions 
are combined on these figures, it can be seen that resolution has no 
discernible effect on 77 as long as the stated resolution standards are 
maintained. 

We consider whether the measured correlations for 77 vs cir¬ 
cular velocity can inform whether energy or momentum are the 
conserved quantities during the blowout. Recall that “energy- 
driven” (energy conserving) winds follow the scaling 77 oc vfl, 
and “momentum-driven/_(momentum conserving) winds follow 
the scaling 77 oc v ^ 1 (Murray et al. 2005). It is expected that low- 
mass galaxies are more affected by energy-driven winds, and high- 
mass galaxies by momentum-driven winds i jMurra y et al.||201lt 
[Hopkins et al.||2012bj ). We find that at high redshift, our low-v c 
halos have a dependence steeper than vfl, while those at high mass 
have a dependence that matches v ^ 1 ■ Nonetheless, a clear transition 
is seen between the two regimes. Stellar feedback on small scales 
in our simulations includes both energy and momentum input from 
several processes in a complex environment. In order to explain the 
dependence of mass-loading factor on halo properties, we there¬ 
fore plan a more detailed study of the nature of wind-driving in our 
simulations in future work. 


5.2 77 at low-z: z < 0.5 

We have elected thus far to show data at z < 0.5, but not to use them 
in our fitting procedure. As we mentioned previously, our sample 
of z < 0.5 non-satellite halos is limited, as it mainly consists of 
the “main” halos from each simulation. Furthermore, the three L*- 
progenitors - ml2i, ml2v, and ml2q - no longer appear to have 
outflows in the CGM that correlate with star formation. We have 
concluded that although, a positive integrated flux of material was 
measured through the 0.25R v i, shell, it is misleading to report this 
as a proper measurement of r /. This is because our methods for 
measuring outflow rates cannot always discriminate between stellar 
feedback-driven winds and other sources of gas motion, such as 
random motion or close-passages of satellites and mergers. 

At epochs z > 0.5, most of the major outflow events coincide 
with high star formation rates seen 50-100 Myr prior. At z < 0.5, 
in the case of the L*-progenitors, star formation rates and outflow 
rates are continuous rather than episodic, and appear to be uncor¬ 
related. Furthermore, using a aw radial velocity threshold to dis¬ 
criminate against random motion generally reduces the outflow rate 
to near-zero values (see Figure]?}, while at earlier times, the reduc¬ 
tion was considerably more modest (see Table [2}. In addition, we 
have visually checked all significant outflow episodes in our three 
L*-progenitors in this redshift interval and found these to originate 
in a close passage of gas rich satellites or the material stripped from 
a passing satellite galaxy. This is in contrast to high-redshift where 
we can connect all significant outflow episodes to the star formation 
bursts. 

Even if we simply include all material with v m d > 0 and make 
a measurement of 77, ml2v and ml2i clearly have 7/ < I, and ml2q 
has 77 < 1 if we exclude the close passage event. A truly discerning 
measurement to only include outflows generated by stellar feed- 


© 0000 RAS, MNRAS 000, 000-000 








14 


Muratov et al. 


4.0 > z > 2.0 



2.0 > z > 0.5 



Figure 7. Top: Total mass of outflows that traversed shells of fixed fraction 
of virial radius between z = 4.0 and z = 2.0 for mlO, mil, ml2v, ml2i, 
and z2h506 respectively. Each value is normalized by the stellar mass of the 
halo at z = 2, which is given in the legend. This quantity generally decreases 
with radius, suggesting that a large fraction of ejected material can either 
stay in the CGM or recycle back to the galaxy, while entrainment is not 
significant. Bottom: The same calculation done for mlO, mil, ml2v, and 
ml2i between z = 2.0 and z = 0.5, normalized by stellar masses at z = 0.5. 
Less material is recycled and retained in the CGM at these intermediate 
redshifts. In the case of mlO, entrainment may play a significant role. 


back would suggest 1 for all L*-progenitors, as can be seen 
from the v cu , = aiD outflow measurement in Figure]?] We find this 
encouraging, as galaxies like the Milky Way that only have low 
to moderate star formation at present day do not appear to typi¬ 
cally launch winds with high mass-loading factor into the CGM 
jVeilleux et al.|2005~]|Rubin et al.|2014| >. This contrasts with the val¬ 
ues of rj that are used for halos of this mass in some large-volume 
simulations which are able form disk galaxies with correct mass, 
but assume that the mass-loading factors remain high at low red- 
shift, as we discuss in the next subsection. 


5.3 Mass-loading as a function of radius within halo 

In order to understand how the flux of material traverses different 
spatial regions of the halo, we provide Figure [7] The top panel of 
this figure shows the total integrated mass of gas that has crossed 
various thresholds in the 4.0 > z > 2.0 interval, normalized by the 
total stellar mass of each halo at z = 2. The values from this plot 
can be used to estimate mass loading at different radii using the 
values of rj at 0.25 R m given in Table]?] 

It is striking that the total gas mass expelled through inner 
region of galactic halos at high redshift is tens or even hundred 
times larger than the total stellar mass accumulated by z = 2. Two 
of the L*-progenitors, ml2v and ml2i, as well as the massive dwarf 
mil, and the LBG-like z2h506 only have about 33% of the total 
material ejected into the CGM (to 0.151Cir) eventually leave the 
halo, while mlO loses nearly all of the outflows. This leads us to 
conclude that at high redshifts, the majority of CGM outflows in 
sufficiently massive galaxies stays within the CGM where a larger 
fraction of this gas can recycle back into the galaxy or contribute to 
the gaseous reservoirs of halos. 

The bottom panel of Figure|7]reveals that the outflow proper¬ 
ties differ at intermediate redshift (2.0 > z > 0.5). In this interval, 
the L*-progenitors and mil now lose 60-70% of material that is 
ejected. This suggests that although their gravitational potentials 
are deepening, these halos are actually more efficient at expelling 
baryons into the IGM. The implications will be discussed further 
in Section [6~2| mlO is particularly peculiar, as it loses significantly 
more mass in the outer regions than the amount ejected from the 
inner regions. The mass of this halo is below the filtering scale in¬ 
duced by the UV background i Thoul & Weinberg] 1996} |Gnedin[ 
|2000||Faucher-Giguere et al.|201 1| , which means that the halo may 

be gradually heated and unbound. Alternatively, this may imply en¬ 
trainment of loosely bound material in the outer region of the CGM 
during the outflow episodes as they propagate outward. We have 
verified that most of the mass loss follows bursts of star formation, 
suggesting that entrainment is the dominant mechanism that pro¬ 
duces this behavior. These calculations do not account for the fact 
that the physical virial radius can be much larger at z = 0.5 than it 
was at z = 2, which means that outflow rates we compute here are 
only valid with respect to the instantaneous position of the virial 
radius but not all of this material actually leaves the growing virial 
radius of a halo. 


5.4 Wind velocity 

As can be seen from Figure]!] our simulations produce winds with a 
broad range of velocities. The detailed kinematic structure of winds 
will be studied in subsequent work, but here we briefly characterize 
the typical wind velocities seen in our simulation. 

We use the same redshift intervals and criteria for inclusion 
as for measurements of r/. At each snapshot, we calculate the 50th 
percentile (median) outflowing (v > Okm/s) radial wind velocity 
from the distribution of flow rates in the 0.25 R W n shell (e.g. left 
panel of Figure]!}. We then average these wind velocities over time, 
weighting the value calculated at each snapshot by the outflow rate. 
This ensures that the average wind velocity is a characteristic of the 
epochs when the most significant outflows are likely to be observed. 
We use this same procedure to also compute the interval-averaged 
95th percentile velocity to give an estimate of some of the fastest 
winds generated by star formation in our simulations. 

We plot the flux-weighted average 50th and 95th percentiles 
in Figure [8] We find strong evidence for a slightly super-linear cor- 
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Figure 8. Flux-weighted average 50th (left) and 95th (right) percentile wind velocity vs v c for all halos in our sample measured at 0.25/? vlr . Same coloring 
convention as Figure[5] Dotted lines show the fits given in Equations |9|and [T()| Both 50th and 95th percentile wind velocities scale slightly superlinearly with 
v c . and are normalized such that typical wind velocities range from ~l-3v c . The shaded region in the left panel shows the range of velocities between best-fit 
lines for 25th and 75th percentile wind velocities. 


relation between wind velocity and v c that is well described by a 
single power law. Correlations between wind velocity and galaxy 
mass are in fact found in observational campaigns that sample a 
sufficiently broad range of galaxy masses ( jMartin et al.|2012) . We 
note that while the 50th percentile velocities are sometimes close 
to the escape velocity of the halo, the 95th percentile velocities are 
considerably faster, and approach 1000 km/s for the most massive 
halos in our sample. Winds with such velocities could be confused 
for winds generated by black holes. It remains to be seen if in mas¬ 
sive compact galaxies that are not present in our sample, stellar 
feedback-driven outflows can be even faster, as observed in some 
extreme cases l |Diamond-Stanic et al.|20 1 2) . 

Again, the three L*-progenitor halos are exceptions at the low- 
z interval, and appear to have significantly slower winds than simi¬ 
lar halos of their mass did at higher redshift. We interpret this as fur¬ 
ther evidence that the outflow rates measured for the L*-progenitors 
at low redshift are not driven by star formation. We find the 50th 
percentile wind velocity is fit by the relation: 

Vwind,50 = 0.85VJ' 1 . (9) 

The fit for 95th percentile wind velocity is fairly similar, dif¬ 
fering mainly in normalization, which is approximately a factor of 
~2 higher: 


v„ind,95 = 1.9V* 1 . (10) 

We have again only used data from the high-z and med-z 
regimes to construct these fits, and consider no redshift evolution, 
as it is not apparent from the figure. We caution that the veloci¬ 
ties discussed here were measured at 0.25R v ir for each halo con¬ 
sidered. Observations of COM winds, whether through quasar ab¬ 
sorption lines or down-the-barrel spectra, cannot easily determine 
the precise distance of outflowing material from the wind’s point 
of origin. This means that a typical observation may include fast- 


moving material within 0.25R v j r that is not considered in our mea¬ 
surements, though it may also include slow-moving material in the 
outer reaches of the halo. In future work, we plan to more closely 
mimic an observational approach and characterize what outflow ve¬ 
locities would be measured at a given SFR and galaxy mass. 

We have verified that the fits are not significantly different 
when we use a v m d > ctid threshold for outflowing gas. This is 
because our flux weighting already filters out material that is tech¬ 
nically outflowing, but too slow to be considered a wind. 


6 DISCUSSION 

6.1 Comparison to subgrid prescriptions 


We present a brief comparison between our results and the subgrid 
prescriptions that have been employed by several other groups (see 
also|Zahid et al.|20l4)|Lu et al.|20 1 5fr. For more details, see Ap¬ 


pendix C We focus on comparing our results to the Illustris project 


l |Vogelsberger et aL| 2013| VI3), as well as the simulations pre- 
sented in |Ford et al.|2014| (F14). Wind velocity measurements from 
the FIRE simulations discussed here are described in Section [54] 
and converted to approximate velocities at the time of wind launch¬ 
ing, using methodology described in Appendix |C| 

The results of our comparison are given in Table [3] There are 
interesting similarities and differences between the wind velocities 
and mass-loading factors found in our work and those from F14 
and VI3. At high redshift, the values of r/ used in VI3 are sys¬ 
tematically higher than ours. The agreement is better for the most 
massive halos in our simulations, but worsens gradually down to 
v c = 60km/s. This is likely because they assume a fixed r) oc v~ 2 
scaling at all masses, while we find r/ ex v7* for v c > 60km/s. The 
launch velocities of winds in V13 are consistent with ours. F14 use 
systematically lower rj at all masses than we do. but launch their 
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Figure 9. Mass fraction of baryons, gas, and stars within R vn relative to the cosmic mean for several runs. The left column shows the three L*-progenitor 
halos (ml2v upper left, ml2i middle left, ml2q lower left). The right column shows mlO (upper right), mil (middle right), and z2h506 (lower right). All 
halos are shown from 3.0 > z > 0 except z2h506, which stops at z = 2. Also shown is the expelled fraction, which is the ratio of M ou t{z), the cumulative mass 
of outflowing material which has traversed 1.0/? v i r , and Mh(z)fbarunh the expected baryonic mass of the halo assuming that it contains an exact proportion of 
the universal baryon budget. This ratio is meant to indicate when the outflows played a significant role in depleting the halo of gas. A low value implies that 
outflows played no significant role in the depletion, while a high value implies that outflows were important. Integration for M ou t(z ) starts at z = 9.0 and is 
only calculated for the most massive progenitor. 


winds at much higher velocities, suggesting that their winds are 
much less likely to ever recycle]^] 

At low redshift, F14 use values of 77 that are in better agree¬ 
ment with our results, while V13 continue to have systematically 
higher values of 77 . For example, VI3 use winds with 77 « 7 at z = 0 


3 The velocities from our simulation quoted in Table pi use an NFW po¬ 
tential to translate CGM velocities discussed in Section l5T4l into launch ve¬ 
locities. We have also tried this exercise using a singular isothermal sphere 
potential, which may be more appropriate for describing massive (L*) halos 
that are baryon dominated in the center at low redshift (Chan et al.|2015) . 
This results in launch velocities up to 30% faster than the NFW results, 
bringing us to closer agreement with F14 for sufficiently massive halos. 


for a Milky Way-mass galaxy (v,„ M ss 200), while our simulations 
find that these halos typically have 77 <C 1. Observations generally 
show that galactic winds are weaker at low redshift (e.g.|Heckman| 
|2001[ [Heckman et al.|20 15 ), which is generally consistent with our 
results, but in tension with V13. F14 uses 77 « 1 for a 10 12 Mq halo. 
As demonstrated in Figure]?] the outflows seen in ml2i at low red¬ 
shift are most likely not related to winds generated by stellar feed¬ 
back, and are generally some combination of random gas motion 
in the halo and close passages of satellites. The values for velocity 
and 77 discussed here should be treated as upper limits. 

Although the winds in our simulations have lower values of 77 
than those of V13, and lower velocities than those of F14, the FIRE 
simulations - like V13 and F14 - nonetheless roughly reproduce 
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the M* — Mi, relation (Hopkins et al.|2Qj~4] >. The key to understand¬ 
ing how this is possible may be the burstiness of star formation in 
the FIRE simulations. Since the consequence of each burst of star 
formation is the dispersal of the ISM, the resultant wind not only 
carries out the gas available for star formation, but also has strong 
dynamical effects on the halo and galaxy. In other words, although 
the galactic wind bulk properties in the FIRE simulations are dif¬ 
ferent than the prescriptions used in V13 and FI4, their detailed 
dynamics (as well as their phase structure) may result in different 
dynamical states for the halo and galaxy, thus limiting the efficiency 
of star formation. 

6.2 Implications for galaxy evolution 

Our analysis has primarily focused on measuring the mass-loading 
factor via outflow rates in the inner regions of the CGM (at 
0.25R v ir). In Section nn we briefly devoted attention to demon¬ 
strating that in the L*-progenitors at high redshift, ml2v and ml2i 
specifically, about a third of material that is ejected into the CGM 
and crosses 0.15R V ir eventually flows out through the virial radius 
of the halo. 

The fact that these numbers are not 100% implies that there 
is a significant amount of material that is initially ejected into the 
CGM. but later able to recycle back into the inner halo and the 
galaxy. Even gas that flows out of the virial radius is not necessarily 
permanently unbound from the halo. This is true by construction, 
as our fiducial choice of v cw = Okm/s is below the local escape 
velocity. In addition, the gravitational potential and virial radius of 
all central halos continue to grow with time, and may re-absorb the 
ejected material. Our results imply that the CGM of massive halos 
hosts a vast reservoir of gas that has been enriched by local high- 
redshift and intermediate-redshift outflows which helps fuel later 
stages of galaxy formation. 

It is interesting to estimate the amount of gas pushed out of 
halos, and to study the evolution of the halo baryon fraction, which 
is a good tracer for the aforementioned reservoir. In Figure [9] we 
show the baryon fraction within R^ over an interval 3.0 > z > 
0 for ml2v, ml2q, ml2i, mil, and mlO, as well as 3.0 > z > 
2.0 for z2h506. We plot the baryon fraction alongside a quantity 
which describes the cumulative mass of outflowing gas which has 
traversed the virial radius since z = 9, relative to the baryon budget 
of the halo if it were to contain the cosmic mean fraction of baryons. 
We call this quantity the expelled fraction. This figure can be used 
to trace the role that galactic winds play in the overall evolution 
of the halo. We note that for many of the halos considered, the 
sum of the baryonic mass fraction within the halo at z = 0 and 
the expelled fraction add up to a number that is near unity. This is 
not necessarily true by construction, as we do not consider whether 
outflows that cross R v i r are actually unbound form the halo, do not 
account for the growth of R v j r through cosmic time, and only follow 
a single massive progenitor. A complete understanding of what sets 
the baryon fraction of a given halo requires us to consider all of 
these processes, as well as any other process that modifies the ratio 
of accretion rates for gas and dark matter. Nonetheless, the expelled 
fraction is a useful diagnostic. 

A low-mass dwarf like mlO already has a depleted baryon 
fraction by z = 3 and it becomes further depleted by a series of pow¬ 
erful outflows between z = 3 and z = 1.5. At lower redshift, outflow 
rates are weaker but infall is only able to replenish a small fraction 
of what has been lost. Infall rates may be diminished due to heating 
of the IGM by the UV background ( jThoul & Weinberg|l996| . The 
sum of the baryon fraction and the expelled fraction are far below 


unity at all times, suggesting that the UV background plays a role 
in preventing accretion continuously, starting at very high redshifts 
(Gnedin(2000}[Faucher-Giguere et al.|2011[ l. The somewhat more 
massive dwarf (mil) can more efficiently accrete the UV heated 
gas, allowing it to have a high baryon fraction while still maintain¬ 
ing powerful outflows until late times (z ~ 0.6), when a merger- 
driven burst expels a huge amount of gas, apparently sealing its 
fate. 

In the L*-progenitor ml2v, star formation is more efficient at 
early times, building up the stellar fraction to -15% of the baryonic 
budget by z ~ 1.5 while maintaining an expelled fraction of un¬ 
der 30%. However, this changes during the interval 1.5 > z > 1.0, 
when the halo undergoes its most prominent burst of star forma¬ 
tion, reaching M* « 40 Mq yr~ 1 . After this powerful burst, as in 
the late-time burst of mil, we see a rapid rise in the expelled frac¬ 
tion. Similar events occur sometime during the med-z interval in 
both ml2q and ml2i. In the case of ml2v, what follows is 7 Gyr 
of evolution where the expelled fraction and baryon fraction stay 
relatively constant. Although accretion does not cease at late times, 
it is slow, and only able to partially repair the damage done by the 
outflows of the intermediate redshift era. This implies that some of 
the ejected material never returns back into the halo that ejected 
it, suggesting that such gas can transport metal polluted gas into 
regions of the IGM that are far from galaxies. 

Following the intermediate-redshift major outflow events, the 
L*-progenitors are left with enough gas to continue star formation, 
but with a somewhat depleted reservoir in the CGM. The three ha¬ 
los exist on a spectrum of their capacity to form stars at late times. 
ml2q has the most efficient star formation at early times, and vi¬ 
olent outflows that did permanent damage to the baryon fraction 
prior to z = 3, followed by more violent activity at intermediate red¬ 
shift. At late times, star formation rates are the lowest of the three 
L*-progenitors, and it is structurally most similar to an elliptical or 
lenticular galaxy (Van de Voort et al., in preparation). ml2v estab¬ 
lishes a stable gaseous disk at late times, but earlier outflows were 
able to deplete the CGM significantly at intermediate redshifts. Fi¬ 
nally, ml2i is able to keep the largest reservoir of gas in the CGM, 
allowing it to have sustained and active star formation at late times. 
We also find that the central galaxy of ml2i has a significant thin 
gaseous and stellar disk at late times. 

The feedback model employed in our simulations starts at 
local scales where the energy and momentum injection rates are 
taken from stellar evolution models and simple approximations to 
feedback processes including supernovae and radiation pressure are 
used to couple the energy and momentum to the ambient gas. Cou¬ 
pling with the local gas properties and galactic conditions then in¬ 
fluences the prevalence, strength, and nature of the emergent phe¬ 
nomenon of global galaxy-scale winds. We have identified that 
variations in halo mass and cosmological redshift of the galaxy at 
the time of wind-launching systematically influence the velocities 
and mass-loading factors of the winds. The underlying mechanisms 
that translate local feedback from star-forming regions to galactic 
winds will be better explored in a future, more focused study. We 
refer readers to |Hopkins et al. | ( j20 14^ for evidence that neither ra¬ 
diation pressure nor supemovae alone are sufficient to reproduce 
the results of the fiducial FIRE feedback model - namely, star for¬ 
mation histories and stellar-mass halo-mass relations. It is evident 
that hydrodynamic processes, gravitational deceleration, as well as 
the multitude of phenomena included in our feedback model affect 
galaxies of different masses and epochs at different efficiencies. 
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Table 3. Comparison of galactic wind mass-loading factors and wind velocities in our simulations to other subgrid models 



halo properties 



wind velocity 


mass-loading 


Name 

z 

M h 

a \/l 

V max 

| Vwind 

Vwind, fit 

Vwind,F 14 

v wind,V 13 

1 v 

Vfit 

Vfu 

7?V13 

mlO 

0.25 

7.1e9 

20 

33 

73-76 

79-99 

170 

80 

45 

54 

28 

270 

mil 

0.25 

1,4e 11 

57 

83 

200-230 

210-270 

410 

210 

7.8 

3.2 

3.5 

42 

ml 2 i 

0.25 

9.8el 1 

120 

180 

340-350 

400-530 

770 

460 

0.31 

1.6 

1.3 

8.9 

mlO 

1.25 

4.3e9 

21 

33 

68-87 

71-95 

180 

84 

110 

89 

26 

270 

mil 

1.25 

l.Oell 

60 

90 

230-310 

200-290 

460 

230 

10 

5.8 

3.1 

36 

ml 2 i 

1.25 

6 . 0 el 1 

100 

150 

350-510 

370-530 

750 

380 

2.3 

3.3 

1.5 

13 

ml 1 

3.0 

2 . 8 el 0 

49 

70 

140-210 

150-230 

410 

180 

14 

14 

4.7 

59 

ml 2 i 

3.0 

1 . 2 ell 

63 

110 

230-330 

250-390 

500 

280 

8.0 

8.7 

2.8 

24 

z2h506 

3.0 

2.3ell 

86 

130 

330-530 

320-490 

670 

330 

6.5 

7.0 

1.7 

17 

z2h506 

2.0 

1 . 2 el 2 

140 

220 

430-600 

500-760 

1100 

560 

6.5 

3.3 

1.0 

6.0 


(1) Name: simulation used 

( 2 ) z: redshift considered 

(3) M/ 7 : halo mass in Mq 

(4) cy\p: radial velocity dispersion within |R V ir as estimated from our simulations (km/s) 

(5) v max \ maximum circular velocity in halo as measured in our simulations (km/s) 

( 6 ) v W i n d'. wind launch velocity from the FIRE simulations, computed by using direct measurement of median and 95th percentile 
wind velocity at 0.25R v j r and re-scaled to launch velocity using the NFW potential (km/s) 

(7) Same as ( 6 ) but median and 95th percentile wind velocity at 0.25/? v ir are computed with our fitting formula (Equations!^ 
and then re-scaled to launch velocities (km/s). ( 8 ) v w i n d } Fu: approximate wind launch velocity computed using|Ford et al.|p014| 
prescriptions (km/s). 

(9) v w i n dyI?,: approximate wind launch velocity computed using|Vogelsberger et al.|J2013| prescriptions (km/s). 

(10) 77 : Direct measurement of mass-loading in FIRE halos 

(11) r)fn : mass-loading factor computed using our fit formula (Equation s[4| and [5j 

(12) rjf 14: approximate mass-loading factor computed using^ 

(13) 77 ^ 13 : approximate mass-loading factor computed using 


Ford et alT](2014) prescriptions 


Vogelsberger et al!||2013| prescriptions. 


7 CONCLUSIONS 

We have explored the relationship between star formation and 
outflows in the CGM in the FIRE hydrodynamic cosmological 
zoom-in simulations. Our key results include: 

(i) We have shown that star formation at high and intermediate 
redshifts (z > 1 ) occurs almost exclusively in episodic bursts. 
The feedback from each burst of star formation drives a "gusty" 
galactic superwind through the CGM. These superwinds are 
sometimes capable of ejecting the entirety of the dense ISM from 
the galaxy and can temporarily suppress gas infall in the inner halo 
(e.g. Figure[2}. However, at lower redshifts we have found that in 
massive halos with M* « 10 12 Mq, star formation becomes less 
bursty and significant outflows largely cease (e.g. Figure]?}. 

(ii) We have quantified 77 - the mass-loading factor - which is 
defined as the ratio of gas mass carried out in galactic winds to 
the mass of stars formed in the galaxy for many of the halos in 
our simulations over specified redshift intervals. At z > 2, we find 
values of 77 as measured in the inner CGM (0.25R v ir) are above 
unity for massive LBG-progenitors, while 77 ~ 10 is typical for 
L*-progenitors, and 77 « 100 are found in low-mass dwarfs. We 
found that 77 varies systematically with halo mass (A/;,), stellar 
mass (M*), circular velocity (v c ), and redshift. The relationship 
between v c and 77 (Figure [5} at low mass is steeper than expected 
for the simplest energy conserving wind model, and at high masses 
it is consistent with the simplest momentum conserving wind 
models (both assume that wind velocities are proportional to v c ). 

(iii) We provide convenient fit functions for 17 and outflow velocity, 
which can be used as a starting point of galactic wind implemen¬ 
tation in large-volume simulations and semi-analytic models (see 


Equations[4|-[T0|and Figures|5]-[8}. 

(iv) Gas throughout the halo is dynamically disturbed following 
particularly powerful outflow events, and star formation can 
be suppressed for significant lengths of time. This regulatory 
mechanism helps to keep halos on the M* — Mi, relation throughout 
cosmic history (see |Hopkins et al.|[2014} . Large-volume simula¬ 
tions which use subgrid models with continuous star formation to 
match the M » — Mi, relation must assume higher values of 77 and/or 
wind velocity than what was found in our work. 

(v) During the intermediate-redshift epoch (2 > z > 1) that 
coincides with the peak of cosmic star formation, the dwarfs 
and L*-progenitors in our sample drive winds that carry material 
beyond the CGM, into the IGM. As cosmic accretion rates begin to 
decrease at these times, many halos can no longer easily replenish 
their supply of baryons. 

(vi) The baryon fraction at low redshift (z < 1) remains low in 
dwarf galaxies, but reaches 40-90% of the cosmic mean for the 
L*-progenitors by z = 0 (Figure|9}. 

(vii) At low redshift (z < 1), dwarfs with a sufficient supply of 
gas can continue to episodically form stars and drive winds with 
77 > 5. At the same time, L*-progenitors no longer drive outflows 
into the CGM through stellar feedback, and stable gaseous disks 
develop. Stars still form at low to moderate rates within these disks, 
but the dense ISM is not dispersed by outflows. This behavior 
demonstrates that our feedback model, which treats early radiative 
feedback from young stars as well as supernovae, is simultaneously 
capable of producing powerful starburst-driven outflows at early 


© 0000 RAS. MNRAS 000, 000-000 


















Gusty, gaseous flows of FIRE 19 


times and the formation of gaseous disks at late times. 

In companion papers, we will address the phase structure of 
the CGM at all times, provide a more detailed history of its 
composition and origin, determine its kinematic structure and 
study observational signatures of galactic winds. 
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APPENDIX A: METHODS FOR CALCULATING 
OUTFLOW RATES 

Al Outflow rate calculation 

Here, we discuss the nuances of the method we used to measure 
the outflow rate, which was introduced in Section [3] In this sec¬ 
tion, we refer to our primary method as the Instantaneous Mass 
Flux method. We describe it in detail, and compare it to another 
method, Inteiface Crossing, which relies on particle tracking. This 
comparison serves as a consistency check for the methods. 

Al.l Method 1: Instantaneous Mass Flux 

If outflow rates were time-steady, and followed a constant velocity 
distribution, this method would robustly capture the outflow rate at 
a given radius. In reality, outflows are time-variable and the veloc¬ 
ity distribution of outflowing particles is unique to each episode. 
This implies outflow rates may not be sufficiently sampled when 
< vsph- Here, At is the time interval between snapshots, and 
vsph is the velocity of an SPH particle in the outflow, and dL is the 
shell length. In this regime, when an outflowing particle that origi¬ 
nated interior to the shell traverses the entirety of the shell between 
snapshots, it is not included as in the rate measurement. On the flip 
side, it is possible to overemphasize the prevalence of fast-moving 
particles if the shell happens to encapsulate a narrow, coherent burst 
at the time of sampling - this transient, high outflow rate will be ex¬ 
trapolated for the entirety of the time interval. Over a large enough 
sample of bursts, these two effects should average out, and average 
measured outflow rates would reflect true values. 

Despite the limitations of this method, we show that it gener¬ 
ally provides accurate estimates of the outflow rate as long as ^ is 
not too far from the typical velocity of SPH particles in outflows. In 
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Figure Al. Cumulative distribution of outflowing particles during a major 
outflow event which crossed the 0.2R v j r threshold since the last timestep 
(black). These particles are counted towards the outflow rate according to 
the Interface Crossing method. We follow these particles backwards in time 
to the previous 5 snapshots from z = 2.04 to z = 2.14, and confirm that 
the majority were within 0. lR v ir two snapshots (—60 Myr) ago. The closest 
halo-centric passage, for each particle within the last 5 snapshots (-150 
Myr) is also shown (dotted line). At the earliest time considered, we see 
that many of the particles now in the outflow were actually accreted form 
outside of 0.1R v ir only to be ejected by stellar feedback. 


our high-redshift regime, the sampling of our simulation snapshot 
interval (~30 Myr) and choice of dL (0. lR v ir, which is 3-5 physical 
kpc for the L* progenitors) would be ideal for sampling particles 
with v ra d between 100 — 200km/s. We find that this is a typical ve¬ 
locity for much of the outflowing material in these halos, though it 
may miss the fastest material. 

Measuring outflow rates with this method gives instantaneous 
results. This is most relevant when contrasted to particle tracking 
methods that apply no velocity cut to particles deemed outflows, 
instead relying only on position differences (e.g. |Ford et al.|20T3) . 
Outflows found through particle tracking methods are not necessar¬ 
ily confined in a narrow range of radii. 

A1.2 Method 2: Interface Crossing 

In this method, particles that are marked as outflows via Equation[2] 
at a given snapshot are traced back exactly one snapshot. We track 
the change in each particle's radial distance from the halo center, 
to determine if it has crossed one of several pre-determined radial 
“interfaces”. We determine the location of the interfaces by using 
the midpoints of the dL = 0.1R v ir shells employed in the Instan¬ 
taneous Mass Flux method. The outflow rate at each interface is 
simply computed by dividing the total mass of all outflowing parti¬ 
cles that passed this interface since the previous snapshot, m sPH , 
by the time elapsed since the previous snapshot. At. Each parti¬ 
cle may cross more than one interface per snapshot, and will be 
counted as part of the outflow rate at every interface it traverses. 
When this method is employed, fast-moving particles are always 
accounted for, regardless of their distance from the halo center. 

The capacity of this method to pick out outflowing material is 
shown in Figure |AT| Here, we see the outflow following a major 
burst of star formation in the central halo of ml2v at z ~ 2 using 


z 



Figure A2. Time evolution diagram for z2h400 from 4.5 > z > 2.0 us¬ 
ing the same format as Figure [5] Red outflow rates and light blue inflow 
rates were computed using the standard Instantaneous Mass Flux method 
(Method 1). We also show the outflow rate at 0.25R v i r and 1.0R v i r as com¬ 
puted by the Interface Crossing Method (Method 2) in dark blue. We see 
that for the most part, the methods produce remarkably similar values for 
outflow rate at both radii. The star formation rate is plotted in black. 


Vcut = cti d- The outflowing particles are picked out by the interface¬ 
crossing method using a threshold at 0.2R v ir- At this epoch, the par¬ 
ticles are mostly contained between 0.2R ylI and 0.4J? v ir, but they all 
resided in the inner regions of the halo and the ISM before begin¬ 
ning their outward trajectory two snapshots (60 Myr) ago. This im¬ 
plies that outflows usually start in the galaxy and ISM, and usually 
do not contain much additional swept-up material from the CGM. 

This method has a few other caveats. As slow-moving and 
fast-moving particles are weighted equally, the velocity cut from 
Equation[2]could influence the measured flux more drastically than 
in the case of the Instantaneous Mass Flux method. In addition, 
because of the specific crossing thresholds, centering errors in the 
halo catalog from snapshot to snapshot can generate spurious sig¬ 
nal, particularly when v cu , = 0 is used. 

A comparison between the outflow rate derived from using In¬ 
stantaneous Mass Flux and Interface Crossing on z2h400 is shown 
in Figure |A2| The similarity between the outflow rates as measured 
by both methods confirms their consistency. The measured values 
of 77 for this halo differ by only 5% between the two methods. 

A2 Alternative velocity cuts 

As mentioned in Section[3] the majority of measurements quoted in 
the text define outflows as all particles with v ra d > v cu t with v cu , = 
0. While using this cut ensures that we account for all gas mass 
that is crossing our specified radii, it may not be the best definition 
for what should count as a galactic wind. Observationally, a wind 
should be distinct from the randomized dispersive motions of gas 
in the halo. For this reason, we provide additional measurements 
of 77 using Vcut = ctid , the 1 -dimensional velocity dispersion of the 
halo. We assume isotropic velocity dispersion such that cr id = 7 ^ 2 . 
We take < 73 d from the halo finder as calculated explicitly using all 
species (gas, stars, and dark matter). This cut provides a simple 
scaling with halo mass that filters out random motion not associated 
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with galactic winds. We note that there is nothing that prohibits gas 
accretion with radial velocity \v ra d\ < o\ d, meaning that the total 
inflow rate may exceed the "infall rate". On the other hand, it is far 
less likely that the halo may launch powerful coherent winds with 
Vrad < cr id, as the local escape velocity in the galactic region will 
prohibit this material from reaching outer regions of the halo. 

In Table]?] we demonstrate that typically, use of this threshold 
reduces measured outflow rates by -25% in L*-progenitors com¬ 
pared to measurements with v cu i = 0. However, the discrepancy 
may is more significant at low redshift: we have argued that r/ is dif¬ 
ficult to measure in the L*-progenitors at z < 0.5 because the out¬ 
flows appear uncorrelated to star formation, and this is supported 
by the low outflow rates of material with v cu1 = aw in Figure]?] 

APPENDIX B: INSTANTANEOUS METHODS FOR 
MEASURING r/ 

B1 Instant rates 

Here, we discuss how we correlate outflow rates and star formation 
rates on the finest possible sampling timescale that is available in 
post-processing of our simulation suite, which is the interval be¬ 
tween snapshots. 

As the star formation rates and outflow rates are highly vari¬ 
able, and are occasionally both measured to be zero, we impose a 
few restrictions on which snapshots we use as data points. First, we 
employ a requirement that the outflow rate must be above zero and 
the specific star formation rate must be above 10~ 10 yr~ 1 . Such cuts 
are necessary to produce a statistically meaningful fit to data, as 
outlying points could differ by several orders of magnitude, and are 
likely the result of mismatched events rather than physical effects. 
We impose one other cut that is somewhat more arbitrary: snap¬ 
shots are only counted if the outflow rate is more than 1 % of the 
maximum measured outflow rate for the given redshift interval in 
the halo. This restriction excludes additional data points produced 
by relatively insignificant outflows set by the halo’s own activity 
pattern. We choose not to impose a symmetric requirement for star 
formation rate, as it excludes too much data when star formation 
is more continuous over an interval. While this could imply that r/ 
is artificially boosted, we find measurements to be consistent with 
other methods. 

Each snapshot’s star formation rate is matched with an outflow 
rate at a time dr later. The timescale dr provides a measurement of 
the delay between the formation of a population of stars, presum¬ 
ably near the halo center, and the time that the bulk of outflowing 
material arrives to a given location in the halo, in this case 0.25f? v ir. 
To understand which value of dr provides the cleanest measure¬ 
ment of 77 , we perform an error analysis on fits for outflow rate 
vs. star formation rate for a range of values of dr (Mout(t + dr) 
vs SFR(t)). We correlate points for each halo separately, over the 
three redshift intervals (high-z 4.0 > z > 2.0, med-z 2.0 > z > 0.5, 
low-z 0.5 > z > 0). 

The results of our correlations for several of the “main” ha¬ 
los in the high-z interval are shown in Figure |B~i~| A range of halo 
masses are represented in this simulation, but we see that dr « 60 
Myr provides the best correlation for all halos shown. Analysis of 
a larger, complete sample including all of the available simulations 
reaffirms that dr « 60 Myr is the mean, median, and mode of "most 
correlated" delay value for halos of all mass in the high-z interval. 
We have verified that this is remarkably close to the expected length 
of time it takes a gas parcel to travel between the halo center and 



d.r (Myr) 



Figure Bl. Top: Reduced x 2 statistic for fits to data points on the M ou ,(t-\- 
dr) vs SFR(t ) relation as a function of dr. The slope of the relation is fixed 
such that Moui{t + dr) oc SFR(t). i.e. linear relation. Each line represents 
one of the “main” halos in various simulations (ml2v, mil, ml2q, and 
z2h600) during the 4.0 > z > 2.0 (high-z) interval. The minimized value of 
X 2 is marked for each halo with a box around the point. In all cases shown 
here, dr = 60Myr is optimal. Bottom: instantaneous outflow rate vs. star 
Formation rate using dr = 60Myr, using the same four halos over the high- 
z interval. While correlation for each halo is apparent, the scatter is still 
high. 


0.25i?vir, if one assumes that the wind travels with constant velocity 
predicted from Equation [9] (62 Myr, 57 Myr, and 54 Myr for mil, 
ml 2 i and z 2 h 600 at z = 3, respectively). Performing the same anal¬ 
ysis for the med-z and low-z intervals, we find best that the average 
preferred value of dr shifts to 90 and 120 Myr, respectively. 

The emergence of this trend with redshift is expected because 
of the growth of physical I? v j r relative to comoving R v i, with red¬ 
shift (we consistently use shells of fixed comoving fraction of R v j r 
when computing outflow rates) for all halos, while the velocity dis¬ 
tribution of outflowing particles during significant bursts does not 
change significantly. The lack of a significant trend for dr as a func¬ 
tion of halo mass suggests that the timescale derived from the ratio 
of R v „ to some characteristic wind speed is constant for all halos, 
and is well represented by dr. We adopt 60 Myr, 90 Myr, and 120 
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Myr as uniform values of dr at 0.25R v ir for the high-z, med-z, and 
low-z regimes, respectively. 

As can be seen from the bottom panel of Figure |BT| the scatter 
about the M OItt (t + dr) vs SFR(t) relation is high even for the op¬ 
timized value of dr. This is indicative of the complex physics that 
goes into wind generation, and shows that each burst of star for¬ 
mation and outflow episode must be studied individually. On top 
of this, we emphasize that our measurements of outflow rate was 
meant to capture material that is flowing into the CGM (which we 
define as 0.25R v ir). In reality, it is possible for outflow rates to be 
high elsewhere in the halo even when it is low at this interface. We 
therefore caution against the interpretation that our optimized val¬ 
ues of dr represent uniquely important timescales on which winds 
should be observed following star formation. 

It is even more difficult to correlate winds in the CGM of the 
outer halo (1 .Of?™) to star-forming events in the galaxy than at the 
inner halo (0.25R v ir). This is because, as we have shown in Section 
|4.1| outflows that reach the virial radius are typically moving more 
slowly than those at 0.25R V ir, and leave the halo at a wider range of 
times. Using a similar \ 2 optimization technique, we find an aver¬ 
age value of dr = 200Myr (during the high-z interval) to be optimal 
for correlating the SFR to outflow rate at 1 .OR™, with no obvious 
trend for halos to have higher or lower values of dr depending on 
mass or v c . 

We used this method to derive average values of 77 for halos 
in our sample (i.e. measure the y-intercept of the best-fit correla¬ 
tion for each halo as shown in the bottom panel of Figure [Erl) , and 
provide measurements in Table [2] 


B2 Episodic Rates 

In this approach, the outflow and star formation rates are integrated 
over an interval sufficient to capture the entire stellar mass formed 
during a single burst of star formation, and the entire mass of the 
gas blown out in the outflow episode. If done properly, this method 
provides nearly the same statistical flexibility as correlating instan¬ 
taneous rates, while also mitigating the stochasticity with using a 
single fixed value of dr. 

However, in order to best identify the episodes over which to 
integrate, we must employ a parametric algorithm that allows for 
flexibility in the duration of the bursts. This episode-finding algo¬ 
rithm must follow the time evolution of the halo (e.g. Figure [3j, 
look for the tell-tale signs that a burst is beginning, integrate the 
star formation for only as long as the burst lasts, then integrate the 
outflow rate for only as long as the outflow is associated with said 
burst of star formation. We constructed such an algorithm, and ver¬ 
ified its validity by ensuring that it was efficient in segregating the 
entire outflow and star formation history into distinct episodes. Pa¬ 
rameters were adjusted to ensure that winds were not misattributed 
to star formation that was unlikely to be physically related to the 
winds and vice versa. Interestingly, we found that the best fits were 
achieved when we mandated that the starting times for integrating 
star formation and outflow rates were separated by timescales com¬ 
parable to the optimal dr measurements with the instant method. 

Using this detection algorithm we find r/ for each episode by 
dividing the total mass of gas expelled per episode by the total mass 
of stars formed. We then again combine the data for all halos in the 
same redshift interval and derive a correlation. We used this method 
to derive average values of 77 over three redshift intervals for halos 
in our sample, and provide the measurements in Table [2] 

The true strength of the instantaneous methods, both instant 
and episodic, could be harnessed by aggregating simulation data 


and statistically analyzing 77 and its dependence on halo properties 
at the epoch of a burst using data points from as many simulations 
as possible. However, we found that our simulation sample was not 
sufficiently large or complete to derive useful results with this ap¬ 
proach. Future work with statistically representative halo popula¬ 
tions, aided by analysis which investigates the origin of the system¬ 
atic scatter, can make optimal use of these measurement techniques 
to thoroughly quantify galactic wind generation. 


APPENDIX C: COMPARISON WITH OTHER WORKS 

Here, we attempt to compare our findings for best-fit values of 77 
and wind velocity to the subgrid models employed in work by other 
groups. We focus on representative stat-of-the-art large-volume 
cosmological simulation projects that did not resolve star-forming 
regions within galaxies, and instead use parametric models to pre¬ 
scribe r? and wind veloc ity. Namely, |Ford et al7||2014| (F14), and 
Vogelsberg er et al. |20 13 1 (V13 Q Both of these works assign wind 
velocities based on halo properties, and then decouple the wind par¬ 
ticles from hydrodynamic interactions to ensure escape from the 
local environment. The following expands on the comparison sum¬ 
marized in Section RxTl and Table[3] 

We emphasize that we cannot precisely reproduce the wind 
prescription employed in these simulations. While each describes 
prescriptions for 77 and wind velocity at the time of launch, the 
results shown in this work have focused on measurements in the 
CGM. Furthermore, the subgrid prescriptions used by those authors 
involve properties such as the local velocity dispersion and stellar 
metallicity, which we did not consider in our fits. Since each group 
used different codes and initial conditions for their simulations, the 
resulting halos are also intrinsically different. As such, we can only 
compare approximations of their models to our results. 

First, we convert our measured wind velocities in the CGM to 
approximate launch velocities near star forming regions. We take 
each outflowing wind SPH particle’s kinetic energy at the time it is 
detected as an outflow at 0.25R v ir and assume that it has conserved 
total energy since being launched at n almc h = min(0.02Rm, lkpc). 
Therefore, we consider the difference in gravitational potential en¬ 
ergy between the two locations. As an approximation for the poten¬ 
tial, we assume that each halo follows an NFW profile with a con¬ 
centration parameter given by Dutton & Maccio (2014). We find 
that launch velocities at high redshift are typically different from 
CGM velocities at 0.25R V ir by a factor of ~2 at the most, but they 
can differ by a factor of -5 or more at low redshift. In addition to 
NFW, we have also considered a singular isothermal sphere (SIS) 
potential, which better describes the inner regions of massive halos 
in the FIRE simulations (e.g. m!2i, fChan et al.|2015| l. We normal¬ 
ize the profile such that the circular velocity is equal for both NFW 
and SIS at r = 0.25R v i r . We have also tried normalizations based on 
the enclosed mass within 0.25R v ii in our simulations. Launch ve¬ 
locities are found to increase by between 10-30%, with the greatest 
increase seen in massive halos that are baryon-dominated at low 
redshift. 

On the other hand, our low-redshift estimate for ml2i is 
also potentially misleading, since the outflow velocity measured 
at 0.25R™ is likely to be random motion rather than a coherent 
star formation-driven galactic wind. As we saw in Section |5.3| 77 
may also differ between 0.25 R v ir and the inner regions of the halo. 
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However, since F14 and VI3 both consider only winds that can es¬ 
cape the innermost regions of halos, and re-couple winds when they 
are outside galaxies, we assume that our measurements of rj in the 
CGM constitute an appropriate analog. 

V13 determine bulk properties of galactic winds using the lo¬ 
cal 1-D dark matter velocity dispersion of the regions where star 
formation takes place. We have no way of reproducing this quan¬ 
tity precisely. However, they suggest that this quantity correlates to 
the halo’s v max according to the work of Okamoto et al. (2010[^ 
We can therefore measure the v max of halos in our simulation and 
approximate o\d.dm • r\ is prescribed based on v„, the wind velocity: 
77 oc v,7 2 . For our comparison, we employ the same parameters used 
in the fiducial model of VI3, which is clearly outlined in that work. 

F14 use the velocity dispersion of the tightly bound "galac¬ 
tic" component of the halo. We approximate the relevant region to 
be i/Jvir, based on the linking length used in their galaxy finding 
algorithm as described in |Oppenheimer & Dave|f2008| l. Wind ve¬ 
locity is set by this velocity dispersion multiplied by several boost 
factors (see |Oppenheimer & Dave|2008[[Ford et al.| 20 l 3| [2014} . 
One of the boost factors depends on metallicity, which we com¬ 
pute using a fit relation for gas-phase metallicity as a function of 
galactic stellar mass for the FIRE simulations l |Ma et al.[2015j ). An¬ 
other is a uniformly sampled random number between 1.05 and 
2, which we simply take to be 1.5. They calculate r/ for winds 
also using the inner velocity dispersion, with a scaling of 77 oc v~ 2 
for v c < 75km/s and rj oc v ^ 1 for v c > 75km/^] The exact feed¬ 
back formula employed in F14 has been provided to us by R. Dave 
(private communication). We refer to the radial velocity dispersion 
within A R v ir as < 77 / 3 , which we measure directly from our simula¬ 
tions, and plug into their formulae. Although velocity dispersion 
fluctuates throughout a given redshift interval, we find our typical 
oja are consistent within 10% of values given by fits fromlHoeftl 

■ m „ „io_ _ __i 1, _ . • 


et al. 


Ford et a 


(20041. Values of 77 are quoted for 10 12 M Q and 10 u Mq halos 


T[(|2014|i, and we roughly reproduce these results. 


We perform the comparison by using three halos in our sam¬ 
ple at z = 3, z = 1-25 and z = 0.25. We consider mlO, mil, ml2i, 
and z2h506. We evaluate z2h506 at both z = 3 and z = 2 for com¬ 
pleteness, as it greatly increases its stellar mass at this interval. We 
measure v max directly from our halo finder to enable the compari¬ 
son. 


5 |Okamoto et al. ( 2010 ) actually suggest that v max is related to the velocity 
dispersion of gas. 

6 It is worth noting that the break in the power law employed in F14 occurs 
at 15km/s which is close to the 60 km/s break we favor in our fit. 
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